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Abstract 

This paper reviews and extends searches for the direct pair production of the scalar super- 
symmetric partners of the top and bottom quarks in proton-proton collisions collected by the 
ATLAS collaboration during tbe LHC Run 1. Most of the analyses use 20 fb”^ of collisions 
at a centre-of-mass energy of = 8 TeV, although in some case an additional 4.7 fb“* of 
collision data at ^[s = 1 TeV are used. New analyses are introduced to improve the sensitiv¬ 
ity to specific regions of the model parameter space. Since no evidence of third-generation 
squarks is found, exclusion limits are derived by combining several analyses and are presen¬ 
ted in both a simplified model framework, assuming simple decay chains, as well as within 
the context of more elaborate phenomenological supersymmetric models. 
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1. Introduction 


In a theory with broken supersymmetry (SUSY) [1-9], the mass scale of the supersymmetric particles is 
undetermined. However, for SUSY to provide a solution to the hierarchy problem [10-13] some of the 
new SUSY particles masses are typically required to be below about one TeV [14,15], hence they could 
be within the reach of the LHC. 

The scalar partners of the right-handed and left-handed chiral components of the top-quark state and 
respectively) are among these particles. In many supersymmetric models, the large Yukawa coupling 
of the top quark to the Higgs sector makes the Higgs boson mass sensitive to the masses of the scalar top 
(referred to as stop in the following) states, such that, to avoid fine tuning, their masses are often required 
to be light. The and components mix to form the mass eigenstates and ( 2 , being defined as fhe 
lighfer of fhe fwo. The scalar superparfner of fhe lefl-handed chiral componenf of fhe bottom quark (b^ 
belongs fo fhe same weak isospin doublel as fhe hence fhey usually share fhe same supersymmefry- 
breaking mass paramefer: a lighf slop can fherefore imply the existence of a light scalar bottom. The 
lightest sbottom mass eigenstate is referred to as 

The ATLAS and CMS collaborations have searched for direct production of stops and sbottoms [16-35] 
using about 4.7 fb“^ of data from the proton-proton collisions produced by the LHC ai ^[s - 1 TeV 
and 20 fb“' at = 8 TeV. These searches have found no evidence of third-generation squark signals, 
leading to exclusion limits in many SUSY models. The aim of this paper is to summarise the sensitivity 
of the ATLAS experiment to R-parity-conserving^ [38-42] models including the direct pair production of 
stops and sbottoms using the full ^[s=^ TeV proton-proton collision dataset collected during Run 1 of the 
LHC.^ The third-generation squarks are assumed to decay to the stable lightest supersymmetric particle 
(LSP) directly or through one or more intermediate stages. The analyses considered are those previously 
published by the ATLAS collaboration on the topic, together with new ones designed to increase the 
sensitivity to scenarios not optimally covered so far. A wide range of SUSY scenarios are studied by 
combining different analyses to improve the global sensitivity. 

The paper is organised as follows: Section 2 briefly reviews fhe expecfed phenomenology of fhird- 
generafion squark producfion and decay; Section 3 reviews fhe general analysis approach followed by 
fhe ATLAS collaboration for SUSY searches; Secfions 4 and 5 presenf fhe exclusion limifs obfained in 
specific models by combining fhe resulfs of several analyses. Two differenl fypes of models have been 
considered: simplified models, where the third-generation squarks are assumed to decay into typically 
one or two different final sfafes, and more complex phenomenological supersymmefric models, where fhe 
slop and sboflom have many allowed decay channels. Conclusions are drawn in Secfion 6. 

For the sake of brevity, the body of the paper provides no details of the ATLAS detector and object 
reconstruction, of the analyses used in the limit derivation, or of how the signal Monte Carlo simulation 
samples were generated. However, a comprehensive set of appendices is provided to supply additional 
information to the interested reader. Appendix A briefly summarises fhe layouf of fhe ATLAS deleclor 
and the general principles used in the reconstruction of electrons, muons, jets, jets containing fj-hadrons 
(fj-jets), and the missing transverse momentum vector (whose magnitude is referred to as L"™®*). 
Appendix B discusses the analyses used to derive the exclusion limits presented in Sections 4 and 5. 


* It is also assumed that the decay of the third-generation squarks is prompt: long-lived and metastable stops/sbottoms are 
discussed elsewhere [36,37]. 

^ The analysis exploiting the measurement of the tt cross section discussed in this paper also uses 4.7 fb“* of proton-proton 
collisions at = 7 TeV. 
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The analyses that have already been published are only briefly reviewed, while those presented for the 
first time in this paper are discussed in detail. Appendix C provides further details of a combination of 
analyses which is performed for the first time in this paper. Finally, Appendix D provides details about 
the generation and simulation of the signal Monte Carlo samples used to derive the limits presented. 


2. Third-generation squark phenomenology 

The cross section for direct stop pair production in proton-proton collisions at = 8 TeV as a function 
of the stop mass as calculated with PROSPINO [43,44] is shown in Figure la. It is calculated to next-to- 
leading order accuracy in the strong coupling constant, adding the resummation of soft gluon emission at 
next-to-leading-logarithmic accuracy (NLO-i-NLL) [45-47]. In this paper, the nominal cross section and 
its uncertainty are taken from an envelope of cross-section predictions using different parton distribution 
function (PDF) sets and factorisation and renormalisation scales described in Ref. [44]. The difference in 
cross section between the sbottom and stop pair production is known to be small [46] , hence the values 
of Figure la are used for both. 



Figure 1: (a) Direct stop pair production cross section at -\/i = 8 TeV as a function of the stop mass. The band 
around the cross section curve illustrates the uncertainty (which is everywhere about 15-20%) on the cross section 
due to scale and PDF variations, (b) Illustration of stop decay modes in the plane spanned by the masses of the 
stop (fj) and the lightest neutralino (Fj), where the latter is assumed to be the lightest supersymmetric particle and 
the only one present among the decay products. The dashed blue lines indicate thresholds separating regions where 
different processes dominate. 

Searches for direct production of stops and sbottoms by the ATLAS collaboration have covered several 
possible final-state topologies. The experimental signatures used to identify these processes depend on 
the masses of the stop or sbottom, on the masses of the other supersymmetric particles they can decay 
into, and on other parameters of the model, such as the stop and sbottom left-right mixing and the mixing 
between the gaugino and higgsino states in the chargino-neutralino sector. 

Assuming that the lightest supersymmetric particle is a stable neutralino (x^), and that no other super- 
symmetric particle plays a significant role in the sbottom decay, the decay chain of the sbottom is simply 
^ bx° (Figure 2a). 
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(e) (f) (g) (h) 

Figure 2: Diagrams of fj and pair production and decays considered as simplified models: (a) b^b^ —» bx\bx^y, (b) 
Vi Z^Jx\tx\\ (c) three-body decay; (d) four-body decay; (e) (f) f/j ^ bx^bx^\ (g) bj}^ tx^txf, 

(h) bj)^ bx^ 2 ^x^- The diagrams do not show “mixed” decays, in which the two pair-produced third-generation 
squarks decay to different final states. 


A significantly more complex phenomenology has to be considered for the stop, depending on its mass 
and on the^j mass. Figure lb shows the three main regions in the plane that are taken into ac¬ 
count. They are identified by different values of = nit^-m^o. In the region where Am^t^,x^^ > 

nit, the favoured decay is —> tx^ (Figure 2b). The region where + mt, < Am^t^,x^^ < mt is char¬ 
acterised by the three-body decay^ (fj ^ Wbx^ through an off-shell top quark, Figure 2c). The region 
where the value of Am{l^,xfj drops below mw + nib, sees the four-body decay bff'x^^, (where / 
and f indicate generic fermions coming from the decay of an off-shell W boson, Figure 2d) competing 
with the flavour-changing decay^ —> cx^ of Figure 2e; the dominant decay depends on the details of 
the supersymmetric model chosen [50]. 

If the third-generation squark decay involves more SUSY particles (other than the j), then additional 
dependencies on SUSY parameters arise. For example, if the lightest chargino (y*) is the next-to-lightest 
supersymmetric particle (NLSP), then the stop tends to have a significant branching ratio for bx\ 
(Figure 2f), or, for the sbottom, b^ —> if kinematically allowed (Figure 2g). The presence of additional 

particles in the decay chain makes the phenomenology depend on their masses. Several possible scen¬ 
arios have been considered, the most common ones being the gauge-universality inspired m^± - 2m^o, 
favoured, for example, in mSUGRA/CMSSM models [51-56]; other interpretations include the case of a 
chargino almost degenerate with the neutralino, a chargino almost degenerate with the squark, or a char¬ 
gino of fixed mass. Another possible decay channel considered for the sbottom is —> bx^ bhx^ 
(Figure 2h), which occurs in scenarios with a large higgsino component of the two lightest neutralinos. 


^ In scenarios that depart from the minimal flavour violation assumption, flavour-changing decays like —> cx° or q —> mF® 

could have a significant branching ratio up to Am (fpF?) ~ 100 GeV [48]. 

The decay —> mF”, in the assumption of minimal flavour violation [49], is further suppressed with respect to —> cF[ by 

corresponding factors of the CKM matrix. 
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(a) (b) (c) 


Figure 3: Diagrams of decays considered as simplified models; (a) (b) t^ki^h', (c) fjfj ^ 

. The diagrams do not show “mixed” decays, in which the two pair-produced third-generation squarks decay 
to different final states. The decay —> yt^ is not an allowed process. 


Despite the lower production cross section and similar final states to fj, the heavier stop state {t^) pair 
production has also been studied: the search for it becomes interesting in scenarios where the detection 
of tj pair production becomes difficult (for example if Am ~ mt). The diagrams of the investigated 

processes are shown in Figure 3. 

Two types of SUSY models are used to interpret the results in terms of exclusion limits. The simplified 
model approach assumes that either a stop or a sbottom pair is produced and that they decay into well- 
defined final states, involving one or two decay channels. Simplified models are used to optimise the 
analyses for a specific final-state topology, rather than the complex (and model-dependent) mixture of 
different topologies that would arise from a SUSY model involving many possible allowed production 
and decay channels. The sensitivity to simplified models is discussed in Section 4. 

More complete phenomenological minimal supersymmetric extensions of the Standard Model (pMSSM 
in the following [57]) are also considered, to assess the performance of the analyses in scenarios where 
the stop and sbottom typically have many allowed decay channels with competing branching ratios. Three 
different sets of pMSSM models are considered, which take into account experimental constraints from 
LHC direct searches, satisfying the Higgs boson mass and dark-matter relic density constraints, or addi¬ 
tional constraints arising from considerations of naturalness. The sensitivity to these models is discussed 
in Section 5. 


3. General discussion of the analysis strategy 

The rich phenomenology of third-generation supersymmetric particles requires several event selections 
to target the wide range of possible topologies. A common analysis strategy and common statistical 
techniques, which are extensively described in Ref. [58], are employed. 

Signal regions (SR) are defined, which target one specific model and SUSY particle mass range. The 
event selection is optimised by relying on the Monte Carlo simulation of both the Standard Model (SM) 
background production processes and the signal itself. The optimisation process aims to maximise the 
expected significance for discovery or exclusion for each of the models considered. 

For each SR, multiple control regions (CR) are defined: they are used to constrain the normalisation of the 
most relevant SM production processes and to validate the MC predictions of the shapes of distributions 
of the kinematic variables used in the analysis. The event selection of the CRs is mutually exclusive with 
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that of the SRs. It is, however, chosen to be as close as possible to that of the signal region while keeping 
the signal contamination small, and such that the event yield is dominated by one specific background 
process. 

A likelihood function is built as the product of Poisson probability functions, describing the observed 
and expected number of events in the control and signal regions. The observed numbers of events in 
the various CRs and SRs are used in a combined profile likelihood fit [59] to determine the expected 
SM background yields for each of the SRs. Systematic uncertainties are treated as nuisance parameters 
in the fit and are constrained with Gaussian functions with standard deviation equal to their value. The 
fit procedure takes into account correlations in the yield predictions between different regions due to 
common background normalisation parameters and systematic uncertainties, as well as contamination 
from SUSY signal events, when a particular model is considered for exclusion. 

The full procedure is validated by comparing the background predictions and the shapes of the distribu¬ 
tions of the key analysis variables from the fit results to those observed in dedicated validation regions 
(VRs), which are defined to be orthogonal to, and kinematically similar, to the signal regions, with low 
potential contamination from signal. 

After successful validation, the observed yields in the signal regions are compared to the prediction. 
The profile likelihood ratio statistic is used first to verify the SM background-only hypothesis, and, if 
no significant excess is observed, to exclude the signal-plus-background hypothesis in specific signal 
models. A signal model is said to be excluded at 95% confidence level (CL) if fhe CL^ [60,61] of fhe 
profile likelihood ratio statistics of the signal-plus-background hypothesis is below 0.05. 

Several publications, targeting specific stop and sbottom final-state topologies, were published by the 
ATLAS collaboration at the end of the proton-proton collision run at = 8 TeV, using a total integrated 
luminosity of about 20 fb“^ Each of these papers defined one or more sets of signal regions optimised for 
different simplified models wifh different mass hierarchies and decay modes for the stop and/or sbottom. 
A few additional signal regions, focusing on regions of the parameter space not well covered by existing 
analyses have been defined since then. All signal regions that are used in this paper are discussed in 
detail in Appendix B, while Table 1 introduces their names and the targeted models. Each analysis is 
identified by a short acronym defined in fhe second column of Table 1. The signal region names of 
previously published analyses are retained, but, to avoid confusion and to ease the bookkeeping, the 
analysis acronym is prepended to their names. Eor example, SRAl from the tOE analysis of Ref. [16], 
which is a search for stop pair production in channels with no leptons in the final state, is referred to as 
tOE-SRAl. 
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Table 1: Summary of the ATLAS analyses and signal regions used in this paper. Each signal region is identified by the acronym of the corresponding analysis 
followed by the original name of the signal region defined either in the published paper or in Appendix B.2. A dash in the signal region name column indicates 
that the analysis does not use the concept of signal region. 
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4. Interpretations in simplified models 


The use of simplified models for analysis optimisation and result interpretation has become more and 
more common in the last years. The attractive feature of this approach is that it focuses on a specific final- 
state topology, rather than on a complex (and often heavily model-dependent) mixture of several different 
topologies: only a few SUSY particles are assumed to be produced in the proton-proton collision - 
often just one type - and only a few decay channels are assumed to be allowed. In the remainder of 
this section, several exclusion limits derived in different supersymmetric simplified models are presented. 
Details about how the MC signal samples used for the limit derivations were produced are available in 
Appendix D. 


4.1. Stop decays with no charginos in the decay chain 

A first series of simplified models is considered. It includes direct stop pair production as the only SUSY 
production process, and assumes that no supersymmetric particle other than the itself and the LSP, 
taken to be the lightest neutralino is involved in the decay. Under this assumption, there is little model 
dependence left in the stop phenomenology, as discussed in Section 2. The stop decay modes are defined 
mainly by the mass separation between the stop and the neutralino, as shown in Figure lb. 

The corresponding diagrams are shown in Figure 2. 

Figure 4 shows the 95% CL exclusion limits obtained in the - m^o plane by the relevant analyses listed 
in Table 1 and discussed in Appendix B, or by their combination. A detailed discussion of which analysis 
is relevant in each range of Am follows. 


Am This kinematic region is characterised by the presence of two competing 

decays: the flavour-violating decay —> cjfj (Figure 2e) and the four-body decay (Fig¬ 

ure 2d). Which one of the two becomes dominant depends on the model details, in particular on the 
mass separation between the stop and the neutralino, and on the amount of flavour violation allowed in 
the model [50]. Several analyses have sensitivity in this region of the - m^o plane. The monojet-like 

signal regions (tc-Ml-3) dominate the sensitivity in the region with Am^t^,xfj ^ mh, regardless of the 
decay of the stop pair, which goes undetected: their selection is based on the presence of an initial-state 
radiation (ISR) jet recoiling against the stop-pair system, which is assumed to be invisible. At larger 
values of signal regions requiring the presence of a c-tagged jet (tc-Cl-2) complement the 

monojet-like signal regions by targeting the —> cx^ decay. Limits on four-body decays can be set using 
signal regions which include low transverse momentum electrons and muons (tlL-bCa_low and WW). 

The limits reported in Figure 4 for these values of Am all assume that the branching ratio of the stop decay 
into either cx^ or tj bff'x^^ is 100%. However, this assumption can be relaxed, and exclusion 
limits derived as a function of the branching ratio of the tj ^ cx^ decay, BR(fj —> cx^), assuming that 
BRffj — > c^°) -I- BR(fj —> bff'x'^) = L Two different scenarios, with Am{t.y,x^^ = 10,80 GeV, are 
considered. The first compressed scenario is characterised by low-px stop decay products, and the set of 
signal regions which have sensitivity is the tc-M, independently of the decay of the stop. In the second 
scenario, the phase space available for the tj decay is larger, and the full set of tc-M, tc-C, tlL-bCa_low, 
tlL-bCa_med and WW-SR selections have different sensitivity, depending on BR(fj —> cx^). 
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Figure 4; Summary of the ATLAS Run 1 searches for direct stop pair production in models where no supersym¬ 
metric particle other than the fj and the is involved in the fj decay. The 95% CL exclusion limits are shown in 
the m,~-m^o mass plane. The dashed and solid lines show the expected and observed limits, respectively, including 
all uncertainties except the theoretical signal cross-section uncertainty (PDF and scale). Four decay modes are 
considered separately with a branching ratio of 100%: fj ^ f^j, where the is mostly fj^, for > m,; 

fj ^ (three-body decay) for mw + mt, < < m,; fj ^ pfj and bff'x^ (four-body decay) for 

Am < mw + mh- The latter two decay modes are superimposed. 


The cross-section limit is derived by combining the analyses discussed above. The SR giving the lowest 
expected exclusion CL^ for each signal model and for each value of BR(f —> cx^) is chosen. Figure 5 
shows the result of these combinations. For Am - 10 GeV, the sensitivity is completely dominated 

by the tc-M signal regions, hence no significant dependence on BR(t —> cx^) is observed. In this case, 
stop masses up to about 250 GeV are excluded. For Am ,x^l'j = 80 GeV, the sensitivity is dominated by 
the tc-C signal regions at high values of BR(f —> cx^). For lower values of BR(r —> cjpj), the “soft-lepton” 
and WW signal regions both become competitive, the latter yielding a higher sensitivity at smaller values 
of the stop mass. The maximum excluded stop mass ranges from about 180 GeV for BR(f ^ - 25% 

to about 270 GeV for BR(f —> cx^) = 100%. 


mw + mt < Am Iii this case, the three-body decay of Figure 2c is dominant. The signal 

regions that are sensitive to this decay are the dedicated signal region defined in the analysis selecting one- 
lepton final states (the tlL-3body) and the combination of several signal regions from the analysis select¬ 
ing two-lepton final states, the t2L. The exclusion limits shown in Figure 4 assume BR(tj —> bWx^) = 1. 
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Figure 5; Upper limits on the stop pair production cross sections for different values of the BRs for the decays 
fj and fj ff’bx^^. Signal points with of 10 GeV (a) and 80 GeV (b) are shown. The limits 

quoted are taken from the best performing, based on expected exclusion CL j, signal regions from the tc-M, tc-C, 
tlL-bCa_low and WW analyses at each mass point. The blue line and corresponding hashed band correspond to 
the mean value and uncertainty on the production cross section of the stop as a function of its mass. The pink 
lines, whose darkness indicate the value of BR(f —» c^j) according to the legend, indicate the observed limit on the 
production cross section. 


The WW signal regions are found to be sensitive to the kinematic region separating the three-body from 
the four-body stop decay region. 


\m In this case, the neutralinos are produced with low pj, and the kinematic properties 

of the signal are similar to those of SM tt production. Exclusion limits in this region were obtained by 
two analyses performing precision SM measurements. The first one is the measurement of the K inclusive 
production cross section cTff. Limits on fj pair production were already set in Ref. [65], which measured 
(Ttt in the different-flavour, opposite-sign channel ep. They were derived assuming a fj decay into an 
on-shell top quark, An extension of the limits into the three-body stop decay is discussed in 

Appendix B. 1 . For a massless neutralino, the analysis excludes stop masses from about 150 GeV to about 
nit. The limit deteriorates for higher neutralino masses, mainly because of the softer b-]e.i spectrum and 
the consequent loss in acceptance. The second analysis considered is that of the top quark spin correlation 
(SC) which considers SM tt production with decays to final states containing two leptons (electrons or 
muons). The shape and normalisation of the distribution of the azimuthal angle between the two leptons 
is sensitive to the spin of the produced particles, hence it allows the analysis to differentiate between stop 
pair and tt production. The limit obtained is shown in the bottom middle (dark orange) of the inset of 
Figure 4. A small region of Am » 180 GeV is excluded with this measurement assuming a small 

neutralino mass. 


Am Iri this kinematic region, the decay tj —> (y-j (see Figure 2b) is dominant. The best 

results in this region are obtained by a statistical combination of the results of the multijet (tOL) and one- 
lepton (tlL) analyses. They both have dedicated signal regions targeting this scenario and the expected 
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sensitivity is comparable for the two analyses. The number of required leptons makes the two signal 
regions mutually exclusive. 



Figure 6: Combined exclusion limits assuming that the stop decays through q ^ fyj with different branching ratios 
X and through fj ^ bxf with branching ratios 1 - x. The limits assume values of x from 0% to 

100% are considered. For each branching ratio, the observed (with solid lines) and expected (with dashed lines) 
limits are shown. 


To maximise the sensitivity to the fj —> decays a statistical combination of the tOL and tlL signal 
regions is performed. The details of the combination are given in Appendix C and the final limit is shown 
in Figure 4 by the largest shaded region (yellow). The expected limit on the stop mass is about 50 GeV 
higher at low m^o than in the individual analyses. The observed limit is increased by roughly the same 
amount and stop masses between 200 GeV and 700 GeV are excluded for small neutralino masses.^ 

A similar combination is performed to target a scenario where the stop can decay as fj —> with 

branching ratio x and as fj —> bxf with branching ratio I - x. Assuming gauge universality, the mass of 
the chargino is set to be twice that of the neutralino. Neutralino masses below 50 GeV are not considered, 
to take into account limits on the lightest chargino mass obtained at LEP [66-70]. The exclusion limits 
are derived for x - 75%, 50%, 25% and 0%.^ Regardless of the branching ratio considered, it is always 
assumed that > rrit + m^o and nr,- > ruh + m^±, such that the two decays t —> tx\ and t —> bx'i are 
both kinematically allowed. A statistical combination, identical to the one described above, is used for 
X = 75%. For smaller values of x, no combined fit is performed, as the sensitivity is dominated by the tlL 
analysis almost everywhere: rather either the tOL or the tlL analysis is used, depending which one gives 
the smaller expected CL^ value. 

Figure 6 shows the result of the combination in the nr,^ - nr^o plane. The limit is improved, with respect 

^ This result holds if the top quark produced in the fj decay has a right-handed chirality. The dependence of the individual 
limits on the top quark chirality is discussed in Refs. [16] and [17], 

^ A value of x = 0% is in fact not achievable in a real supersymmetric model. Nevertheless, this value has been considered as 
the limiting case of a simplified model. 
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to the individual analyses, by about 50 GeV for m^o = 50 GeV and ;c = 75%. For other a: values, the tlL 

analysis is used on the full plane, with the exception of the point at the highest stop mass for m-o = 50 

■^1 

GeV at X - 50% and 25%. Stop masses below 500 GeV are excluded for m,-o < 160 GeV for any value 

r ^1 

of X. 


4.2. Stop decays with a chargino in the decay chain 

In the pMSSM, unless the higgsino-gaugino mass parameters are related by Mi ju, M 2 , the mass dif¬ 
ference between the lightest neutralino and the lightest chargino cannot be too large. The mass hierarchy 
m^o < m^± < is, hence, well motivated, leading to the decay chain shown in Figure 2f. 

If additional particles beside the stop and the lightest neutralino take part in the stop decay, the stop phe¬ 
nomenology quickly becomes complex. Even if the chargino is the only other relevant SUSY particle, the 
stop phenomenology depends on the chargino mass, on the stop left-right mixing, and on the composition 
of the neutralino and chargino in terms of bino, wino and higgsino states. 

Figure 7 shows the exclusion limits obtained by the analyses listed in Table 1 and discussed in Appendix B 
if a branching ratio of 100% for t bx^ is assumed. The exclusion limits are presented in a number 
of nit-m^o planes, each characterised by a different hypothesis on the chargino mass. For all scenarios 
considered, the chargino is assumed to decay as x\ where the (*) indicates a possibly virtual 

W boson. 


= 5,20 GeV. This scenario assumes that the difference in mass between the lightest char¬ 
gino and the neutralino is small (Figure 7a), which is a rather common feature of models where, for 
example, the ESP has a large wino or higgsino component. Two hypotheses have been considered, with 
Am(y*,^j) = 5 GeV and Am(xf,x^) = 20 GeV. Eor both, the complete decay chain is 7j —> bx^ —> 
bff'x^, where the transverse momenta of the fermions / and /' depend on Am(x^,x^^) and on the stop 
mass, given the dependency on the chargino boost. If Amixf’X^) - 5 GeV, the fermions have momenta 
too low to be efficiently reconstructed. The observed final state then consists of two Z>-jets and This 
final state is the direct target of the bOE signal regions. Eor Am(x^,x^) = 20 GeV, the signal efficiencies 
of the bOE signal regions decrease because of the lepton and jet veto applied. The tlE signal regions with 
soft leptons, instead, gain in sensitivity, profiting from the higher transverse momentum of the fermions 
from the off-shell W decay produced in the chargino decay. 


m^± = 106,150 GeV. This scenario (Eigure 7b) assumes a fixed chargino mass. The SR yielding the 

lowest expected exclusion CE^ for this scenario depends on the value of Amfjf'*,^®). Eor Am(x'^,x^) < 20 
GeV, the bOE signal regions provide the best sensitivity; for larger values of Am(x\,x\)’ the tlE and t2E 
signal regions provide better sensitivity because of the same mechanism as in the Am(Y*,^j) = 5,20 
GeV scenario above. The exclusion extends up to about 600 GeV for small values of Am(x^,x\)- A 
region of the parameter space with my- up to about 260 GeV and between 100 GeV and m^± is not 
yet excluded. 
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Figure 7; Summary of the ATLAS Run 1 searches for direct stop pair production in models where the decay mode 
?! ^ bxi with ^ is assumed with a branching ratio of 100%. Various hypotheses on the fj, and 

mass hierarchy are used. Exclusion limits at 95% CL are shown in the fj - x^ mass plane. The dashed and 
solid lines show the expected and observed limits, respectively, including all uncertainties except the theoretical 
signal cross-section uncertainty (PDF and scale). Wherever not superseded by any a/s = 8 TeV analysis, results 
obtained by analyses using 4.7 fb ' of proton-proton collision data taken at Vs = 7 TeV are also shown, with the 
corresponding reference. The four plots correspond to interpretations of (a) the bOL and tlL soft-lepton analyses 
in two scenarios (Am(x^,x^^) = 5 GeV in light green and Am(xf,x'^) - 20 GeV in dark green), for a total of four 
limits; (b) the bOL, tlL and t2L analyses in scenarios with a fixed chargino mass = 106 GeV (dark green) and 
= 150 GeV (light green); (c) the tlL and t2L analyses in scenarios with - 2ot^o; (d) interpretations of the 

tlL, t2L and WW analyses in senarios with Am (fpVf) - 10 GeV. 
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/Mp± = 2 /m~o. Inspired by gauge-universality considerations, the third scenario (Figure 7c) is charac- 
'^1 

terised by a relatively large The t2L signal regions dominate the sensitivity for ~ m^±. 

The sensitivity of the dedicated tlL-bC is dominant in a large region of the plane, and determines the 
exclusion reach for moderate to large values of Am 


= 10 GeV. The fourth scenario (Figure 7d) assumes a rather compressed tj - x^ spectrum. 
The region at low m,-^ and large m^o is characterised by low mass separations between all particles in¬ 
volved, and it is best covered by the tlL-bCc_diag, the tlL soft lepton, and the WW signal regions. At 
larger values of the stop mass, the leptons emitted in the decay have larger px. and the t2L signal 
regions provide the best sensitivity. 



Figure 8: Exclusion limits assuming that the stop decays through b+x'f b + +x^ with branching ratio 

of 100% assuming a fixed stop mass of m,- = 300 GeV. The region below the purple line and above the blue line, 
indicated by a light shading, is excluded. 


nii = 300 GeV. The final scenario considered is one where the stop mass is fixed at 300 GeV, and 
1 

the exclusion limits are expressed in the m^±-m^o plane. In the case of the compressed scenario, corres¬ 
ponding to a small mass difference Am(Yp^j), the fermions from the decay can escape detection 
and only the two (^-jets and F™®® would be identified in the final state. Thus, the bOL signal regions are 
expected to have a large sensitivity in this case, while for larger values of Am(Y*,V?)> the lepton can be 
observed, yielding a final-state signature investigated by the tlL soft-lepton signal region. A combination 
of the bOL and tlL signal regions is performed by choosing, for each point of the plane, the SR giving 
the lowest CL^ for expected exclusion. The result, reported in Figure 8, shows that a large portion of the 
plane is excluded, with the exception of a region where the mass separations between the , the x\ and 
the ^j are small. 

Summarising, in the simplified models with fj ^ bx'i —> bW^*^x^^, stop masses up to 450-600 GeV are 
generally excluded. Scenarios where Am(fj,^j) is small are particularly difficult to exclude and in these 
compressed scenarios, stop masses as low as 200 GeV are still allowed (Figure 7b). A small unexcluded 
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area is also left for a small region around = (180,100,50) GeV (Figure 7c), where the 

sensitivity of the analyses is poor because the signal kinematics are similar to SM tt production. 


4.3. Limits on pair production of 

Although the pair production of 7^ has a cross section larger than that of and although the decay 
patterns of the two particles can be similar, it can be convenient to search for the latter in regions where 
the sensitivity to the former is limited. This is the case, for example, in the region where Am (tpA?) ~ 
of Figure 4, where the separation of 7j pair production from SM top quark pair production is difficult. The 
t2tlZ and t2tlh analyses are designed to detect pair production in this region of the - m^o plane, 
followed by the decays —> t^Z and > t^h. The Higgs boson h is assumed to have a mass of 125 GeV 

and SM branching ratios. 

The exclusion limits were first derived in a scenario in which the pair-produced ^2 decays either through 
^2 —> with a branching ratio of 100% (Figure 3a), or through ^2 —> ht^ (again with a branching ratio 
of 100%; Figure 3b). In both cases, the 7j is assumed to decay through 7j —> and its mass is set to 
be 180 GeV above that of the neutralino (assumed to be the LSP), which is the region not excluded in 
Figure 4. The final state contains two top quarks, two neutralinos, and either two Z or two h bosons. 



Figure 9: Exclusion limits at 95% CL in the scenario where pair production is assumed, followed by the decay 
t 2 Zt^ (blue) or i^h (red) and then by tx\ with a branching ratio of 100%, as a function of the and 
F'j mass. The fj mass is determined by the relation m,- - = 180 GeV. The dashed lines indicate the expected 

limit and the solid lines indicate the observed limit. 


Figure 9 shows the exclusion limits for the t2tlh and the t2tlZ analyses. In both cases, a limit on is 
set at about 600 GeV for a massless neutralino. In the case of a decay through a Higgs boson, the limit 
covers neutralino masses lower than in the case of the decay through a Z boson. 

The assumption on the branching ratio of the ^2 has also been relaxed, and limits have been derived 
assuming that the three decays —> Zt^, Tj —> ht.^ and t 2 —> (Y? (Figure 3c) are the only possible ones. 
The limits are shown in Figure 10 as a function of the three BRs, for different combinations of the ^2 and 
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ATLAS Vi = 8 TeV, 20 fb ' 

tg-tg production, tg Zt,, ht^, 
Am(t^, X°) = 180 GsV, t, -> tx° 


I Observed t2t1Z 
I Observed t2t1h 
I Observed tO/tl L comb. 


Expected t2t1Z 
Expected t2t1h 
Expected t0/t1 L comb. 



Figure 10; Exclusion limits as a function of the branching ratio for t 2 t^h, t^Z and tx\- The blue, 
red and green limit refers to the t2tlZ, t2tlh and combination of tOL and tlL analyses respectively. The limits are 
given for three different values of the and masses. 


masses. Three analyses have been considered: the t2tlZ, t2tlh and the comhination of the tOL and tlL 
discussed in Section 4.1.^ The three analyses have complementary sensitivities. Together, they exclude 
pair production with a mass of 350 GeV and 500 GeV for - 20 GeV. A non-excluded region appears 
for - 500 GeV if larger masses are considered. 

4.4. Sbottom decays 

Under the assumption that no supersymmetric particle takes part in the shottom decay apart from the 
lightest neutralino, the shottom decays as bx^ with a branching ratio of 100% (Figure 2a). The 

final state arising from sbottom pair production hence contains two fj-jets and The bOL signal 

regions were explicitly optimised to be sensitive to this scenario. In case of a mass degeneracy between 
the sbottom and the neutralino, the general consideration that the monojet-like tc-M selection is almost 
insensitive to the details of the decay of the produced particles still holds: the tc-M signal regions offer 
fhe besf sensifivify for scenarios where ~ m^o. 

Figure 11 shows fhe limifs of fhe fc and bOL analyses on fhe - m^o plane. The monojef-like (fc-M) 
SRs exclude models up fo a value of ~ m^o ~ 280 GeV. Sboffom masses are excluded up fo abouf 
600 GeV for neufralino masses below abouf 256 GeV 


’ For the combination of the tOL and tlL analyses, the limits extracted for the f, —» tx° decay with branching ratio of 100% 
have simply been rescaled by appropriate factors depending on the branching ratio of fj —> tx° considered here. 
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Figure 11; Observed (solid lines) and expected (dashed lines) 95% CL limits on sbottom pair production where the 
sbottom is assumed to decay as bx\ with a branching ratio of 100%. The purple lines refer to the limit of the 
tc analysis, while the blue lines refer to the bOL analysis. 


If other supersymmetric particles enter into the decay chain, then multiple decay channels would be 
allowed. Similarly to the stop, the case in which other neutralinos or charginos have a mass below the 
sbottom is well motivated. The branching ratios of the sbottom to the dilferent decay channels depend on 
the supersymmetric particle mass hierarchy, on the mixing of the left-right components of the sbottom, 
and on the composition of the charginos and neutralinos in terms of bino, wino, and higgsino states. 

An exclusion limit is derived under the assumption that the sbottom decays with a branching ratio of 100% 
into —> tx\ (Figure 2g). The chargino is assumed to decay through x^ —> with a branching 

ratio of 100%. The final state is a complex one, and offers many handles for background rejecfion: if 
pofenfially confains up fo fen jefs, fwo f^-jefs, and up fo four leptons. The limifs of Figure 12a, shown in 
fhe - m^Q plane, were obfained by using fhe fhree-lepfon signal regions SS3L, eifher fixing fhe mass 
of fhe neufralino fo m^o = 60 GeV or by making fhe assumpfion fhaf m^± = 2m In fhe fwo scenarios 
considered, sboffom masses up to abouf 440 GeV are excluded, wifh a mild depenciency on fhe neufralino 
mass. 

The lasf case considered is one where fhe pair-produced sboffoms decay fhrough bx^, followed by 
fhe decay of x^^ info a x^ and a SM-like Higgs boson h (Figure 2h). The final sfafe confains up to six 
fj-jets, four of which are produced by the two Higgs bosons decays. Since multiple f7-jets are present in 
the final sfafe, fhe fhree-f?-jefs signal regions (g3b) are used fo place limifs in fhis model. 

The limif, derived as a funcfion of and assuming a fixed neufralino mass of x^ = 60 GeV, is 
shown in Figure 12b. Sboffom masses befween abouf 300 and 650 GeV are excluded forjf^ masses above 
250 GeV. 
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m~ [GeV] m- [GeV] 


(a) (b) 

Figure 12: Exclusion limits at 95% CL for a scenario where sbottoms are pair produced and decay as (a) 
with a BR of 100% or (b) bx^ with a BR of 100%. The signal regions used in (a) are the SS3L, and two 
different models are considered: a fixed neutralino mass of 60 GeV (in purple) or = 2m^o (in blue). The limits 
are shown in the plane. The signal regions used in (b) are the g3b-SR-0j. A fixed neutralino mass of 60 

GeV is assumed, and the limit is shown in the -rn^o plane. 


5. Interpretations in pMSSM models 

The interpretation of the results in simplified models is useful to assess the sensitivity of each signal 
region to a specific topology. However, this approach fails to test signal regions on the complexity of 
the stop and sbottom phenomenology that appears in a realistic SUSY model. To this extent, the signal 
regions are used to derive exclusion limits in the context of specific pMSSM models. 

The pMSSM [57] is obtained from the more general MSSM by making assumptions based on experi¬ 
mental results: 

- No new source of CP violation beyond the Standard Model. New sources of CP violation are 
constrained by experimental limits on the electron and neutron electric dipole moments. 

- No flavour-changing neutral currents. This is implemented by requiring that the matrices for the 
sfermion masses and trilinear couplings are diagonal. 

- First- and second-generation universality. The soft-SUSY-breaking mass parameters and the trilin¬ 
ear couplings for the first and second generation are assumed to be the same based on experimental 
data from, e.g., the neutral kaon system [71]. 

With the above assumptions, and with the choice of a neutralino as the LSP, the pMSSM adds 19 free 
parameters on top of those of the SM. The complete set of pMSSM parameters is shown in Table 2. A full 
assessment of the ATLAS sensitivity to a scan of the 19-parameters space has been performed in Ref. [72]. 
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Table 2: Description of the 19 additional parameters of the pMSSM model with a neutralino LSR 


Parameter 

Description 

wThr, wrjyj, niqii , ifigR, 

First- and second-generation common mass parameters 

fHfRi 

Third-generation mass parameters 

My, M2, M3 

Gaugino mass parameters 


Trilinear couplings 

fi,MA 

Higgs/higgsino mass parameters 

tanyS 

Ratio of vacuum expectation values of the two Higgs doublets 


Here, a set of additional hypotheses are made, to focus on the sensitivity to a specific, well-motivated set 
of models with enhanced third generation squark production: 

- The common masses of the first- and second-generation squarks have been set to a multi-TeV scale, 
making these quarks irrelevant for the processes studied at the energies investigated in this paper. 
This choice is motivated by the absence of any signal from squark or gluino production in dedicated 
SUSY searches performed by the ATLAS [62, 63,73-76] and CMS [29, 34,77-82] collaborations. 

- All slepton mass parameters have been set to the same scale as the first- and second-generation 
squarks. This choice has no specific experimenfal or fheorefical mofivafion, and should be regarded 
as an assumpfion. 

- A decoupling limif wifh Ma - 3 TeV and large fanyS values (fan/3 >15) has been assumed. This is 
partially mofivafed by resulfs of fhe LHC searches for higher mass Higgs boson sfafes [83, 84]. 

- For lan/3 » 1, fhe Higgs boson mass depends heavily on fhe producf of fhe slop-mass paramefers 

fhe mixing belween fhe lefl- and righl-handed slates X, = At - ji! fan/3 [85]. 
The slop seclor is Iherefore complefely fixed, given fhe Higgs boson mass, fhe value of X, and one 
of fhe Iwo slop mass paramefers^. 

- The Irilinear couplings A}, in fhe sbolfom seclor are found lo have limiled impacl on fhe phenomen¬ 
ology, and are Iherefore sel lo zero. 

- The gluino mass paramefer M 3 is sel such fo evade LHC consfrainls on gluino-pair production. 

These assumpfions reduce fhe number of addilional free parameters of fhe model lo fhe mass paramel- 
ers of fhe eleclroweak seclor ip., M \, M 2 ) and Iwo of fhe fhree Ihird-generalion squark mass paramefers 
{mqL^,mtR,m^£). All fhe assumpfions made eifher have a solid experimenfal basis, or are intended fo 
simplify fhe inlerprelalion in terms of direct production of stops and sbottoms (as, for example, the as¬ 
sumption on the slepton mass parameters). 

Three types of models have been chosen, that, by implementing in different ways constraints arising 
from naturalness arguments and the dark-matter relic density measurement, further reduce the number 
of parameters to be scanned over. They are described below, and summarised in Table 3 together with 
additional information on the most relevant production and decay channels. 


In particular, a minimum value of Ms ~ 800 GeV is allowed if the maximal mixing condition XJMs - V6 is realised. 
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Table 3: Details of parameters scanned in the three pMSSM models used for interpretations. The settings of additional parameters that are hxed for each model 
are also given together with the main production and decay channels targeted. 
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Naturalness-inspired pMSSM: The model is inspired by naturalness criteria, which require a value of 
lu in the range of a few hundred GeV, favour stop masses below one TeV, place weak constraints on the 
gluino mass and give no constraints on the mass of other SUSY particles [ 86 ]. The exclusion limits are 
determined as a function of the higgsino mass parameter /r and the left-handed squark mass parameter 
The parameter is scanned in the range 350 GeV < < 900 GeV. The parameter is 

scanned in the range 100 GeV < lu < rriqi^ - 150 GeV, where the lower bound is determined by limits 
on the chargino mass arising from LEP [66-70]. The right-handed stop mass parameter mfR and the stop 
mixing parameter Xt are determined by choosing the maximal mixing scenario XdMs - V 6 and by the 
requirement of having a Higgs boson mass of about 125 GeV. The other squark and slepton masses, as 
well as the bino mass parameter Mi, are set to 3 TeV. The wino mass parameter M 2 is set such that 
M 2 = Sju. The gluino mass parameter M 3 is set to 1.7 TeV. 

With this choice of the model parameters, the spectrum is characterised by two light neutralinos 

and one chargino (uf )> with masses of the order of fj., a light with a mass of the order of and 
a light tj with mass of the order of niqiji up to niqi^, ~ 700 GeV (the constraint on Mj does not allow 
the mass of fj to increase beyond about 650 GeV). The production processes considered are direct pair 
production of and with similar masses. Because of the abundance of light higgsino states, many 
different decays can occur. 


Well-tempered neutralino pMSSM: The models are designed to loosely satisfy dark-matter thermal- 
relic density constraints (0.09 < < 0.15, where h is the Hubble constant), while keeping fine tuning 

(defined as in Ref. [87]) fo less fhan 1%. The exclusion limifs are determined as a function of Mi and 
^qLi, or Ml and nifR, wifh ji ~ -M\ in bofh cases fo safisfy fhe dark-mailer conslrainls Ihrough Ihe 
presence of well-lempered neulralinos [ 88 ]. The conslrainls on Ihe Higgs boson mass are satisfied in a 
way similar fo Ihe naluralness-inspired pMSSM model above. All olher parameters are Ihe same as in Ihe 
naluralness-inspired pMSSM model. These models tend lo have Ihree neulralinos and Iwo charginos wilh 
masses lower lhan fj or giving rise lo a diverse phenomenology. 


/r/Z-enriched pMSSM: These models are defined such lhal Higgs and Z bosons are produced abund- 
anlly in Ihe SUSY particles’ decay chains. The assumption of Mi = 100 GeV ensures Ihe presence of a 
bino-like neulralino LSP, while M 3 = 2.5 TeV ensures lhal direcl gluino production is highly suppressed 
compared lo Ihird-generalion squark production. Two sels of models have been defined: in Ihe firsl one, 
/r and Ihe righl-handed sbollom mass parameter are scanned while keeping M 2 = /r, niqi 2 - 2 TeV, 

rtijR = 1.6 TeV; in Ihe second one, ju and niqLT, are scanned while keeping M 2 = 1 TeV, = 3 TeV, 
ifitR = 2 TeV. The former is dominated by sbollom pair production, while bolh sbollom and slop pair 
production are relevanl for Ihe laller. Slop mixing parameters are chosen wilh maximal mixing lo satisfy 
Higgs boson mass conslrainls. In Ihese models, Ihe decays of Ihe Ihird generation squarks into Ihe heavier 
neulralino slates and^®) are followed by decays to Ihe lighlesl neulralino wilh Ihe emission of a Z or 
a h boson. Typically the ^2 (T 3 ) decays into a Z boson 30% (85%) of the times, and into a Higgs boson 
70% (15%) of the times. The subsequent decays of the Higgs boson into Ir-quark pairs (happening with 
the same branching ratio as in the Standard Model) lead to final slates rich in Ir-jels. 


Exclusion limils for Ihese pMSSM models are determined by combining many of Ihe SRs defined for Ihe 
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searches discussed in this paper (tOL, tlL, tb,^ t2tlZ, g3b, tc). For each set of parameters the individual 
95% CL expected limit is evaluated. The combined exclusion contour is determined by choosing, for 
each model point, the signal region having the smallest expected CL^ value of the test statistic for the 
signal-plus-background hypothesis. 


b, t pair production: Naturai pMSSM 



Figure 13: Expected and observed 95% CL exclusion limits for the naturalness-inspired set of pMSSM models from 
the combination tOL, tlL and tb analyses using the signal region yielding the smallest CLj value for the signal- 
plus-background hypothesis. The dashed black line indicates the expected limit, and the yellow band indicates 
the +lcr uncertainties, which include all uncertainties except the theoretical uncertainties in the signal. The red 
solid line indicates the observed limit, and the red dotted lines indicate the sensitivity to ±lcr variations of the 
signal theoretical uncertainties. The dashed and dotted grey lines indicate a constant value of the stop and sbottom 
masses, while the dashed light-blue line indicates a constant value of the neutralino mass. 

Figure 13 shows the exclusion limit for the naturalness-inspired set of pMSSM models based on the tOL, 
tlL and tb analyses. The tOL and tlL analyses have a similar expected sensitivity. These SRs were 
optimised assuming a 100% BR for or —> bx"^, while for these pMSSM models, the stop 

decays to fj —> tjfj, fj —> bx\ and fj —> bx^ with similar branching ratios (and the sbottom to both 
b^ tx^)- The tb signal regions, discussed in detail in Appendix B.2.3, are designed to be 
sensitive to final states containing a top quark, a f^-quark and missing transverse momentum and address 
such mixed-decay scenarios by requiring a lower jet multiplicity. 

The signal regions that dominate the sensitivity are the tb, tOL-SRCl and tlL-bCd_bulk at low values 
of niqL^, and tb, tOL-SRAl, tOL-SRA2 and tlL-tNbC_mix at intermediate and high values of m^/, 3 . The 
excluded region for models with ~ 900 GeV and yu ~ 150 GeV is due to the saturation of m;- at high 
mgL 3 values: to satisfy the Higgs boson mass constraint requires Ms ~ 800 GeV, hence niy- at nigts ~ 900 
GeV is smaller than that at ~ 800 GeV. The large fluctuations of the observed limit with respect 
to the expected one are due to transitions between different signal regions providing the best expected 
exclusion in different regions of the plane. 


^ The tb signal region, discussed in detail in Appendix B.2.3, implement a one-lepton selection, designed to be sensitive to final 
states containing a top quark, a i-quark and It complements the selections of the tOL and flL signal regions targeting 

final states. 
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b, t pair production: pMSSM with weil tempered LSP t pair production: pMSSM with weii tempered LSP 



Figure 14: Expected and observed 95% CL exclusion limits for the pMSSM model with well-tempered neutralinos 
as a function of M\ and (a) niqui or (b) mj^. The limit of (a) is obtained as the combination of the tOL, tlL, tb and 
SS3L analyses, while the tOL analysis is used for (b). The signal region yielding the smallest CL^ value for the 
signal-plus-background hypothesis is used for each point. The dashed black line indicates the expected limit, and 
the yellow band indicates the +lcr uncertainties, which include all uncertainties except the theoretical uncertainties 
in the signal. The red solid line indicates the observed limit, and the red dotted lines indicate the sensitivity to 
+ 1(T variations of the signal theoretical. The dashed and dotted grey lines indicate a constant value of the stop and 
sbottom masses, while the dashed light-blue line indicates a constant value of the neutralino mass. 


Figures 14a and 14b show the exclusion limit obtained for the set of pMSSM models with well-tempered 
neutralinos as a function of niqui and rtijn, respectively. In both cases, the exclusion is largely dominated 
by the tOL analysis. For Figure 14a, the signal region dominating the sensitivity at low is tOL-SRCl, 
while at higher values tOL-SRAl and tOL-SRA2 dominate the sensitivity. The drop in sensitivity 
at niqi-^ - 410 GeV, M\ = 260 GeV is due to the opening of the and transition, 

kinematically suppressed for smaller values of the difference mqL 3 - M^. Such decays introduce more 
intermediate states in the decay, effectively reducing the transverse momenta of the final state objects. 
The large fluctuations of the observed limit are again due to transitions between different signal regions. 
For Figure 14b, the sensitivity is entirely dominated by the various tOL-SRC. The difference in sensitivity 
between these two scenarios is due to the presence of both a stop and a sbottom for small rriqi^, while 
only a stop is present for low values of mg?. 

Finally, Figures 15a and 15b show the exclusion limit obtained for the set of /i/Z-enriched pMSSM 
models. These models yield large ^-jet multiplicities to the final state through direct sbottom decays, top- 
quark decays and ^2 The exclusion is dominated by the tOL and g3b analyses for Figure 15a 

and by and the tOL analysis for Figure 15b. 

More informations about the limits obtained, including the SLHA files for the points mentioned in Table 3, 
can be found in Refs. [89] and [90]. 
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b, t pair production: pMSSM with M^=100 GeV, M2=1000 GeV 
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Figure 15: Expected and observed 95% CL exclusion limits for the set of h/Z-enriched pMSSM models as a function 
of fi and (a) m ^/,3 and (b) . The limit of (a) is obtained as the combination of the tOL, g3b, t2tlZ and SS3L 

analyses, while the tOL, t2tlZ and tb analysis are used for (b). The signal region yielding the smallest CL, value 
for the signal-plus-background hypothesis is used for each point. The dashed black line indicates the expected 
limit, and the yellow band indicates the +lcr uncertainties, which include all uncertainties except the theoretical 
uncertainties in the signal. The red solid line indicates the observed limit, and the red dotted lines indicate the 
sensitivity to +lcr variations of the signal theoretical. The dashed and dotted grey lines indicate a constant value of 
the stop and sbottom masses, while the dashed light-blue line indicates a constant value of the neutralino mass. 


6. Conclusions 


The search programme of the ATLAS collaboration for the direct pair production of stops and sbottoms 
is summarised and extended by new analyses targeting scenarios not optimally covered by previously 
published searches. The paper is based on 20 fb“^ of proton-proton collisions collected at the LHC by 
ATLAS in 2012 at a centre-of-mass energy ^fs - ^ TeV. Exclusion limits in the context of simplified 
models are presented. In general, stop and sbottom masses up to several hundred GeV are excluded, 
although the exclusion limits significantly weaken in the presence of compressed SUSY mass spectra 
or multiple allowed decay chains. Three classes of pMSSM models, based on general arguments of 
Higgs boson mass naturalness and compatibility with the observed dark-matter relic density have also 
been studied and exclusion limits have been set. Large regions of the considered parameter space are 
excluded. 
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A. The ATLAS detector and object reconstruction 


The ATLAS detector [91] consists of inner tracking devices surrounded by a superconducting solenoid, 
electromagnetic and hadronic calorimeters and a muon spectrometer immersed in a toroidal magnetic 
field. The inner detector (ID), in combination with a superconducting solenoid magnet with a central 
field of 2 T, provides precision fracking and momenfum measuremenfs of charged parficles in a pseu- 
dorapidify^° range |? 7 | < 2.5. The ID consisfs of a silicon pixel defector, a silicon microsfrip defecfor and 
a sfraw lube Iracker (|? 7 | < 2 . 0 ) lhal also provides Iransifion radiation measuremenfs for eleclron idenlific- 
alion. A high-granularily eleclromagnelic calorimeler sysfem, wilh acceplance covering I 77 I < 3.2, uses 
liquid argon (LAr) as fhe active medium. A scinfillalor-file calorimeler provides hadronic coverage for 
\r]\ < 1.7. The end-cap and forward regions, spanning 1.5 < \r]\ < 4.9, are inslrumenfed with LAr elec¬ 
tromagnetic and hadronic calorimeters. The muon spectrometer has separate trigger and high-precision 
tracking chambers which provide trigger coverage for \r}\ < 2.4 and muon identification and momentum 
measurements for \t]\ < 2.7. 

The data sample used in this analysis was taken during the period from March to December 2012 with 
the LHC operating at a pp centre-of-mass energy of = 8 TeV.'^ Following requirements based on 
beam, detector conditions and data quality, the complete dataset corresponds to an integrated luminosity 
of 20.3 fb“\ with an associated uncertainty of 2.8%. The uncertainty is derived following the same 
methodology as that detailed in Ref. [92]. Events used in the analyses presented in this paper were 
selected using the ATLAS three-level trigger following different chains based on the signatures being 
considered. A common set of cleaning cuts, aimed at rejecting events heavily contaminated by non¬ 
collision backgrounds, or events containing badly measured or fake jets is applied to all analyses. 

The experimental signature of third-generation supersymmetric particles includes the production of f^-jets 
in association with missing transverse momentum and possibly additional jets and charged leptons. Dif¬ 
ferent signatures are investigated in this paper to gain sensitivity to a variety of possible topologies arising 
from the production and decay of stops and sbottoms. Different event selections share common definitions 
of the final reconsfrucfed objecfs, which are defailed in the remainder of this Appendix. Analysis-specific 
departures from those definitions are detailed for each case in Appendix B or in the specific analysis 
paper. 

The reconsfrucfed primary verfex [93] is required to be consisfenf with the luminous region and to have 
at least five associafed fracks with pj > 400 MeV; when more than one such vertex is found, the vertex 
with the largest summed Pj of the associated tracks is chosen. 

Jets are constructed from three-dimensional clusters of noise-suppressed calorimeter cells [94] using the 
anti-kf algorithm [95-97] with a distance parameter R = 0.4 and calibrated with a local cluster weighting 
algorithm [98]. An area-dependent correction is applied for energy from additional proton-proton colli¬ 
sions based on an estimate of the pileup activity in a given event using the method proposed in Ref. [99]. 
Jets are calibrated as discussed in Ref. [100] and required to have pt > 20 GeV and I 77 I < 4.5. Events 


ATLAS uses a right-handed system with its origin at the nominal interaction point (IP) in the centre of the detector and the 
z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the ^-axis points upward. 
Cylindrical coordinates (r, are used in the transverse plane, being the azimuthal angle around the beam pipe. The 
pseudorapidity is defined in terms of the polar angle 6 as rj = -lntan(0/2). The distance AR in the r]-(p space is defined as 
AR= ^l(Aijf + (A0)2. 

** The limits derived using a measurement of the tt production cross section discussed in Section 4.1 and Appendix B.l also 
uses 4.7 fb“' of pp collisions data collected at a centre-of-mass energy of ^7 = 7 TeV. 
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containing jets arising from detector noise, cosmic-ray muons, or other non-collision sources are removed 
from consideration [100]. 

Jets arising from a ft-quark fragmentation and within the acceptance of the inner detector (It/I < 2.5) 
are identified with an algorithm that exploits both the track impact parameters and secondary vertex 
information [101]; this algorithm is based on a neural network using the output weights of the IP3D, 
JetFitter-i-IP3D, and SVl algorithms (defined in Refs. [102,103]). A lower cut on the output of the neural 
network defines the (r-tagged jets. Three different working points are used, with a nominal efficiency 
of 60%, 70% and 80% as evaluated on simulated top quark pair production events. The corresponding 
rejection factors against jets originating from light (c) quarks are 25 (3), 135 (5) and 600 (8). 

Electrons are reconstructed from energy clusters in the electromagnetic calorimeter matched to a track in 
the inner detector [104] and are required to have \r]\ < 2.47. Several criteria, including calorimeter shower 
shape, quality of the match between the track and the cluster, and the amount of transition radiation 
emitted in the TRT detector, are used to define three selections with decreasing efficiency and increasing 
purity, named respectively ’loose’, ’medium’ and ’tight’ [104]. These three electron selections are used 
throughout this paper in the definitions of various signal and control regions. Muons, which are identified 
either as a combined track in the muon spectrometer and inner detector systems, or as an inner detector 
track matched with a muon spectrometer track segment [105,106], are required to have \r}\ < 2.4. 

Electrons and muons (generically referred to by the symbol £) are usually required to have transverse 
momentum px > 10 GeV. Eor specific scenarios with compressed mass spectra, low-px leptons are 
expected and the px threshold is lowered to 6 GeV for muons and to 7 GeV for electrons. 

The missing transverse momentum p™*® (with magnitude E™®®) is the negative vector sum of the px 
measured in the clusters of calorimeter cells, which are calibrated according to their associated recon¬ 
structed object (e.g. jets and electrons), and the px of the muons. Calorimeter cells not associated with 
any reconstructed object are also used in the calculation of p™*®. The missing transverse momentum from 
the tracking system (denoted by with magnitude is computed from the vector sum 

of the reconstructed inner detector tracks with px > 500 MeV and \r]\ < 2.5, associated with the primary 
vertex in the event. 


B. Analyses used in the paper 


Several signal regions are used in this paper, either standalone or in combination with others, to derive ex¬ 
clusion limits in the many models considered. This Appendix provides a review of the already published 
analyses and a more extended documentation of the signal regions not previously published. 


B.l. Review of already published signal regions 

The discussion of analyses that have already been published is reduced to a summary for the sake of 
brevity. Table 1 provides a reference to the papers where full details of the signal, control and validation 
region selections, together with the strategies adopted for the estimation of the background processes are 
found. 


28 


Multijet final states (tOL) The analysis is designed to be sensitive to final states arising from all- 
hadronic decays of directly pair-produced stops [16]. Two sets of signal regions were optimised to max¬ 
imise the sensitivity to topologies arising from tj —> decays, assumed to happen with a branching 

ratio of one. The first set of signal regions, named tOL-SRA, assumes that both top quark hadronic decays 
can be fully resolved by indentifying the six final-state jets. The SM background (dominated by tt and 
Z+ heavy flavour (HF) Jefs production) is rejected based on fhe presence of fwo hadronic sysfems con- 
sisfenf wifh fop quarks and large The second sef of signal regions, named fOL-SRB fargefs a similar 
scenario, buf aims af topologies where fhe fop quarks have a large boosf, and some of fhe decay producfs 
are merged info a single jef. The even! selection is designed fo selecf final sfafes wifh a maximum of five 
R - 0.4 anfi-k( jefs, fo be mufually exclusive wifh fOL-SRA, and relies on fhe presence of R = 0.8 and 
R = 1.2 anfi-ki jefs confaining fhe hadronic decay producfs of fhe fwo fop quarks. The jef masses, fhe 
fransverse mass of fhe R™*® and fhe nearesf and ofher variables are used fo discriminafe againsf fhe 
dominanf SM tt, Z + HF jefs and W + HF jefs producfion background processes. 

Finally, a fhird sef of signal regions, named fOL-SRC, is designed fo increase fhe analysis sensifivify fo 
fhe decay fj ^ bx^. The presence of fhe intermediate chargino sfafe fends fo decrease fhe jef mulfipli- 
cify: fhese signal regions require five anfi-k( jefs wifh R = 0.4, and base fhe signal selection on a sef of 
fransverse mass variables aimed af rejecfing fhe dominanf SM tt producfion process. 


One-lepton final states (tlL) The large number of signal regions defined in fhis analysis sfems from 
fhe variefy and complexify of fhe possible stop final sfafes considered [17]. All signal regions are char- 
acferised by fhe presence of one lepfon, a second-leplon vefo, a minimum of fwo jefs and large R™®®. A 
firsl sef of four signal regions (llL-fN) were opfimised assuming a branching rafio of 100% for fhe decay 
t —> tXy These signal regions aim af having sensifivify fo differenf Am(t,x^), in parficular flL-fN_diag 
fargefs scenarios wifh small Ani(fj,^j) and makes use of fhe shape informafion of fhe R™®* and nij disfri- 
bufions.^^ The flL-fN_boosf SR fargefs models wifh fhe largesf Am{t,x\), where fhe fop quark produced 
by fhe sfop decay has a large boosf and large-R jefs are used fo reconsfrucf fhe fop quark decays. 

The decay t Xi^^ infroduces additional degrees of freedom in fhe decay. The final-sfafe kinemafics is 
largely driven by fhe mass separafion befween fhe slop and fhe chargino Am(t,^*), and by lhaf belween 
fhe chargino and fhe neufralino Am{x^,x\). Several signal regions, idenfitied by fhe prefix llL-bC were 
designed and opfimised depending on fhe mass hierarchy and, consequenlly, on fhe differenf kinemafics 
of fhe lepton and f7-jels. 

The four signal regions llL-bCa_low, llL-bCa_med, llL-bCb_medl and llL-bCb_high largel small val¬ 
ues of Am(xf,x^) and have fhe common fealure of making use of a dedicafed soff-lepfon selection: muons 
and elecfrons are idenfitied down fo a pj Ihreshold of 6 GeV and 7 GeV, respectively, requiring a special 
frealmenf for fhe esfimafe of possible background processes arising from lepfon misidenliticafion. They 
are collecfively referred fo as “soff-leplon” signal regions. Bofh flL-bCa signal regions require a hard 
ISR jef fo boosf fhe stop pair system and produce a sizeable R™®®. The flL-bCb fargefs large values of 
Am{t^,X^) and exploifs fhe presence of fwo relatively hard ti-jefs in fhe even!. 


The transverse mass mj of the lepton with transverse momentum pi and the missing transverse momentum vector with 
magnitude is defined as 

mT = ^2(|;fT|Rr'-/fT-pr*) (1) 

and it is extensively used in one-lepton final states to reject SM background processes containing a W boson decaying lepton- 
ically. 
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The signal region tlL-bCc_diag targets a mass hierarchy complementary to that of the tlL-bCb. The 
small value of Am(t,xj) gives rise to soft Z?-jets that go undetected, hence Zj-tagged jets are vetoed for this 
region. 

Topologies arising from scenarios where both Am{t,Xi) and Am(xf,x^) are sizeable are targeted by the 
three tlL-bCd regions: they all require four jets in the final state, are characterised by different f^-jet 
multiplicities, and apply different selections on the mj and amT 2 ^^ variables. A veto on additional 
isolated tracks and r lepton candidates identified wifh loose criferia helps fo suppress fhe dominanf SM 
background from dilepfonic tt decays. 

The lasf fwo signal regions listed in Table 1, flL-3body and flL-lNbC_mix, were opfimised for fwo 
addifional possible scenarios. If Am(t,x^) < nit and fhe mass hierarchy or fhe model paramefers suppress 
fhe decay fhrough a chargino, fhen fhe dominanf sfop decay is f —> bWx^^, fhrough an off-shell lop 
quark (Ihree-body decay). The dedicaled signal region relies on fhe shape informalion from fhe m-j and 
amT :2 variable disfribufions. Finally, flL-lNbC_mix is designed fo recover sensifivify in scenarios where 
fhe slop is assumed fo decay wifh similar probabilities to tjf-j and bx\'- the selection aims to reject the 
dominant dileptonic ft background by making use of the topness [109] variable. 


Two-lepton final states (t2L) If the SUSY mass hierarchy forbids the presence of sleptons in the stop 
decay chain, final slales conlaining fwo leplons {e or jj) and a large amounl of would arise from slop 
pair producfion. The main background is given by SM processes conlaining fwo W bosons in fhe final 
slale (mainly tt and WW) [18]. To discriminafe fhe slop signal from fhe SM background, fhe slransverse 
mass variable mT 2 [107, 108] is used. The slransverse mass, computed using the two leptons as visible 
particles and the missing transverse momentum vector, exhibits a kinematical end-point at for most 
SM processes. Because of the presence of additional due to the LSP, the end-point for a SUSY 
signal can be at larger values, depending on the mass separation between the particles involved in the 
decay. The analysis is optimised assuming t^ X\i> with BR=100% and Am(x^,x^^) > mw, but it is also 
sensitive to the three-body decay mode of the stop. To derive exclusion limits, five signal regions (f2L) 
have been defined, requiring differenl jef mulfiplicifies and differenl mT 2 Ihresholds. A selecfion requiring 
fwo fj-jefs and based on mj 2 compufed using Ihem as visible parficles is sensilive fo fhe chargino decay 
mode wifh Am{t^,x^) > nit. Finally, a mullivariale discriminanf is builf which largels fhe fj —> decay 
mode. 


Final states from compressed stop decays (tc) If the difference in mass between the stop and the 
neutralino is smaller than the W boson mass, then the only possible decay channels are f —> jfjC or 
t W*b, where the decay products of the off-shell IT* would, in general, be soft. This analysis [19] 
has defined two sets of signal regions, both optimised for the t x^c decay. A common preselection 
requires the presence of a high-pj jet, large F™®® and applies a lepton veto. The first set of signal regions 
named tc-M, targets scenarios with the stop mass almost degenerate with the neutralino mass, and applies 
a selection that exploits a monojet-like signature arising from the presence of an ISR jet. Three different 
signal regions have been designed, characterised by increasing thresholds on the leading jet pj and 
The second set of signal regions, named tc-C, targets less compressed scenarios, and exploits the presence 


The asymmetric stransverse mass variable is a variant of the stransverse mass variable [107,108] defined to efficiently reject 
dileptonic tt decays. It assumes that the undetected particle is the W boson for the branch with the lost lepton and the neutrino 
is the missing particle for the branch with the observed charged lepton. For the dileptonic tt events, amt 2 is bounded from 
above by the top quark mass, whereas new physics can exceed this bound. 
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of jets originating from the fragmentation of c-quarks in the final state. A dedicated c-tagging algorithm 
was used to reject the dominant SM background processes arising mostly from tt and Z ^ vv (produced 
in association with heavy-flavour jets) production. As in the case of the tc-M signal regions, different 
thresholds on the leading jet pj and on are used to identify a looser and a tighter tc-C region. 

Final states with a Z boson (t2tlZ) A Z boson can be emitted in the decay of t 2 t^Z, producing 
final sfafes wifh large lepfon mulfiplicifies. If can be useful fo look for (rafher fhan tj) producfion if, for 
example, fhe mass of fj is very close fo fhe sum of fhe fop quark and neufralino masses, which would lead 
fo tj pair production final sfafes difficulf fo disfinguish from SM tt producfion. Models are invesfigafed 
wifh = 180 GeV wifh fhe decay fj —> (y-j. The final sfafe would confain, beyond fhe Z boson, 

several jefs arising from fhe fj decay. Similar final sfafes can be obfained in GMSB models where fhe Z 
boson is emitted in fhe ^ j —> GZ decay if fhe gravifino G is fhe LSP and fhe neufralino fhe NLSP. 

This analysis [20] defines five differenf signal regions divided info fwo sefs. The firsl sef, named l2llZ- 
SR2, requires fwo same-flavour lepfons whose invarianf mass is consisfenf wifh fhaf of a Z boson, mz, and 
af leasf one (^-fagged jef. The fhree signal regions are characferised by fhe differenf selection fhresholds 
applied fo fhe F™*®, fo fhe transverse momentum of the dilepton system pj {£() and to the jet multiplicity. 
The second set of signal regions, named t2tlZ-SR3, requires three leptons, two of which must form an 
opposite-sign same-flavour pair whose invariant mass is consistent with mz- Both signal regions require 
at least five jets, among which at least one has to be (^-tagged. The two signal regions are characterised 
by the different selection thresholds applied to px (^0 and to the leading lepton px. 


Final states with two b-jets and (bOL) This signature arises naturally from the sbottom decay 
Moreover, one expects the same final state from t^ followed by —> ff'x^ in the 

limit of small ^m{x^,x\)- This analysis [21] defines two sets of signal regions, bOL-SRA and bOL-SRB, 
targeting scenarios with large and small squark-neutralino mass separations, respectively. 

The event selection of bOL-SRA requires large exactly two f?-jets and vetoes the presence of addi¬ 
tional jets; the rejection of the SM ft production background is carried out by making use of the contrans- 
verse mass [110] of the two Z^-jets. Its distributions shows a kinematical end-point at about 135 GeV for 
tt production, while extending to higher values for the signal. 

A selection relying on the presence of an ISR jet is instead needed if the third-generation squark mass 
is almost degenerate with that of the neutralino. This is the purpose of bOL-SRB, which selects a hard, 
non-fj-tagged leading jet recoiling against the squark pair system. The selection includes the requirement 
of two f7-tagged jets, a veto on additional hadronic activity, and the presence of large 


Final states with three A-jets (g3b) This analysis [62] is designed to search for gluino-mediated sbot¬ 
tom and stop production in events with no leptons or one lepton (electron or muon) in the final state. 
However, it was found to have sensitivity for direct production followed by ^ x\b —> x^^hb, where 
h is the SM Higgs boson with mass m/, = 125 GeV, and also sensitivity to some of the pMSSM mod¬ 
els considered in this paper. Such final states are characterised by a large multiplicity of f7-jets both in 
gg —> tj^tt and gg ^ b^b^bb where there are up to four (j-jets in the final state. 


31 


Three sets of signal regions have been designed to target different mass hierarchies of the gluino-mediated 
sbottom and stop production models. All signal regions have at least four jets with pj > 30 GeV, three 
identified b-jets, large and a large rtisff, defined as the scalar sum of the pj of the jets and 

Strongly produced final states with two same sign or three leptons (SS3L) Final states containing 
many leptons or same-sign (SS) leptons can arise from the pair production of gluinos and squarks, when 
the produced particles decay to the LSP through multiple intermediate stages, or when several top quarks 
appear as part of the decay chain. The analysis was developed for the gluino-mediated stop production 
process gg —> followed by fj —> which can yield final states containing up to four leptons, 

including SS pairs. Similar final states arise from the sbottom decay —> tx^, which are studied in this 
paper. 

This analysis [63] concentrates on final states containing either three leptons or a SS lepton pair produced 
in association with many jets. Five signal regions (identified by the prefix SS3L) are defined, which are 
characterised by different light- and heavy-flavour jet multiplicities, high selection thresholds on and 
/Weff, and different thresholds on the transverse mass of the lepton with the highest transverse momentum 
and the 

Spin correlation in tt production events (SC) If the mass of the is such that ~ m^o +mt, the final- 
state kinematics are similar to that of Standard Model tt production. One possible approach is to derive 
exclusion limits on the stop mass by performing SM precision measurements. This analysis has measured 
the azimuthal angle difference between the two leptons arising from the dileptonic tt decay [64]. The 
events are required to contain, beside the two leptons, at least two additional jets, one of which is required 
to be f7-tagged. In events containing two leptons of the same flavour, the Z production background is 
suppressed by applying a selection on the dilepton invariant mass. The distribution of the azimuthal angle 
between the two leptons is sensitive to the spin correlations of the tt system: it is hence used to extract 
limits on possible contaminations from direct scalar top production events. 


ti production cross section (xsec) The measurement of the tt production cross section using events 
containing two different-flavour leptons ep and Z?-tagged jets is used in Ref. [65] to extract limits on the 
direct pair production of fj with mass close to the top quark. The assumed decay is fj —> (y-j. 

The tt production cross section cr^- is obtained by using the equations 


Ml - La,Jee^,2et{^ - Ci,eb) + (2) 

N2 = Latjee^Cbel + ivf ® (3) 

where N\ and N 2 are the number of events with two different flavour leptons having exactly one or two 
f?-tagged jets, respectively, L is the integrated luminosity, egfi the efficiency for a tt event to pass the 
lepton selection, is the probability of having a fj-jet within acceptance and for it to be tagged, Cb is a 
correlation coefficient which is close to unity, and and are the number of events with one or 
two f?-tagged jets from SM events different from tt production. The values of (Ttt and et are extracted 
from the data by solving the two simultaneous equations (2) and (3), avoiding the need to estimate et 
from simulation. 
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Stop-pair production events with nif^ > nit + have similar and b-]&i kinematics to SM tt produc¬ 
tion events, so the fitted value of et in a combined sample is compatible with that from tt production 
events alone, and the fitted cross section corresponds closely to the sum of tt and stop-pair production 
cross sections. Limits on stop pair production are extracted by calculating 95% CL limits on the stop pair 
production signal strength /r (defined as fhe rafio of fhe obfained stop cross secfion to fhe fheorefical pre¬ 
diction) based on fhe comparison of fhe measured cross secfion wifh fhaf predicted for SM tt production 
evenfs alone. A 95% CL signal sfrengfh smaller fhan unify for a given signal poinf implies ifs exclusion. 

This inferprefafion, which made use of collision dafa wifh bofh ^/s = 7 and 8 TeV, is exfended here to fhe 
fhree-body decay fj —> Wbx^. The main difference wifh respecf fo fhe scenario considered in Ref. [65] 
is fhaf fhe fhree-body decay fends fo yield f^-jels wifh lower pj, leading fo a tiffed et for fhe combined 
sample which is differenl from fhaf expecfed for tt evenfs alone. The limifs obfained are summarised in 
Figure 16 for a neufralino mass of 1 GeV. A 95% CL limif fhaf excludes slop masses below 175 GeV 
is obfained. The figure also shows fhe effecl on fhe limif of a “sneaky top squark” scenario [111]: fhe 
presence of a t^ wifh mass similar fo fhaf of fhe fop quark could bias fhe measuremenf of fhe lop-quark 
mass ifself. The bias in fhe lop-mass measuremenf inlroduced by fhe existence of a fj wifh mass = 170 
GeV depends on fhe analysis technique and channel, and was evalualed fo be al mosl 1 GeV for fhe Iwo- 
and Ihree-dimensional lemplafe lechniques used in fhe ATLAS fop mass measuremenf in fhe leplon-tjels 
channel [112]. The effecl of a potential bias of 1 GeV and 2.5 GeV on fhe lop-mass measuremenf was 
sludied by recalculaling fhe observed 95% CL limif on p when reducing fhe predicfed SM tt production 
cross secfion from fhe baseline value of nit = 172.5 ±1.0 GeV to Ihose obfained for top mass cenlral 
values of 173.5 and 175 GeV. The corresponding limif on fhe stop mass is reduced by aboul 5 and 15 
GeV, respectively. 

The dependence of fhe exclusion limifs on fhe neufralino mass was sludied and found to be imporfanl: fhe 
effecl of an increasing neufralino mass is fo decrease fhe pj of fhe ti-jefs, and hence fo lower fhe value of 
€b for fhe slop pair production signal. For a neufralino mass of 30 GeV, only a small range of stop masses 
around 150 GeV is excluded. 

The sensilivily of fhe tt cross-section measuremenf to fj pair production assuming a branching rafio of 
100% info 7j ^ bx^, followed hy x\ wifh m^± = 2m^o was also investigated. The presence of 

fhe intermediate chargino slate fends to lower fhe pj of fhe leptons and of fhe ti-jels signiticanlly, hence 
decreasing bofh and et- No exclusion limif can be derived for Ihis scenario. 

Summarising, fhe limifs on slop pair producfion obfained in Ref. [65] have been extended by considering 
fhe stop fhree-body decay. Stop masses belween 150 GeV and nit can be excluded for a neufralino mass 
of 1 GeV. The exclusion holds provided fhaf any bias in fhe fop-quark mass measuremenf by a nearby 
stop is nol signiticanl. Sludies indicate fhaf fhis pofenlial bias would affecl fhe limif on fhe slop mass by 
less fhan 5 GeV. 


B.2. Description of the new signal regions 

New analyses were developed fo largel lopologies and regions of fhe SUSY parameler space nol well 
covered by previously published signal regions. They are idenlitied Ihroughoul Ibis paper and in Table 1 
wifh fhe acronym WW, lb and l2llh. Their confribulion to fhe exclusion limifs derived bofh in simplified 
and pMSSM models is ouflined in Secfions 4 and 5 respecfively. In fhis Appendix, furfher defails aboul 
Ihese analyses are provided for fhe inferesfed reader. Additional informafions aboul seleclion efficiencies, 
sensifivilies of fhe differenl signal regions and individual limif plols, please refer fo Refs. [89] and [90]. 
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Figure 16: Expected and observed 95% CL limits on the signal strength ji (defined as the ratio of the obtained stop 
cross section to the theoretical prediction) for the production of fj pairs as a function of m,-. The stop is assumed to 
decay as tx\ or through its three-body decay depending on its mass. The neutralino is assumed to have a mass 
of 1 GeV. The black dotted line shows the expected limit with +lcr uncertainty band shaded in yellow, taking into 
account all uncertainties except the theoretical cross-section uncertainties on the signal. The red solid line shows 
the observed limit, with dotted lines indicating the changes as the nominal signal cross section is scaled up and 
down by its theoretical uncertainty. The short blue and purple dashed lines indicate how the observed limits with 
the signal cross section reduced by one standard deviation of its theoretical uncertainty for m,- < m, when the top 
quark mass is assumed instead to be 173.5 + 1.0 and 175.0 ±1.0 GeV. 


B.2.1. Final states with two leptons at intermediate values of /mt 2 (WW) 

The measurement of the production cross section of nonresonant WW pairs in the two-lepton channel 
at the LHC [113-115] has given rise to theoretical speculations [116-118] which interpret the possible 
excess as due to the production of a light stop. The mass hierarchy favoured by these speculations includes 
a fj with mass around 200 GeV, a xj degenerate with it, and m^± - m^o of a few tens of GeV: possible 
hadronic decay products of the fj —> bx^ transition would have low pj and would allow the events 
to survive the tight jet-veto selections applied in the SM cross-section measurement. Dedicated signal 
regions, defined by requiring two different-flavour opposite-sign leptons in the final states, are designed 
to have maximum sensitivity to such scenarios. The approach is also sensitive to scenarios where the stop 
decays predominantly through the three-body bWx\ or four-body fj —> b€vx^^ decay. 

MC simulated events are used to model the signal and to describe all backgrounds that produce two 
prompt leptons from W,Z or h decay. For processes whose predicted yield in the signal regions is small, 
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or whose topology resembles very closely that of the signal, making it hard to define a proper control 
region, the background estimate is fully based on MC simulation. For tt, Z + jets and WW production 
processes, which are the dominant backgrounds, the acceptance of the signal regions selection is estimated 
with MC simulation, while the normalisation is estimated in dedicated control regions. The MC samples 
used are the same as in Ref. [18]. 

The identification criteria for electrons, muons and jets follow the strategy defined in Appendix A: 
baseline electrons, which are used in the estimation of the fake-lepton background, are selected by ap¬ 
plying the “medium” identification criteria. Signal electrons are identified using the “tight” criteria, and 
they are further required to be isolated. Signal muons correspond to baseline muons with an additional 
calorimeter- and track-based isolation requirement applied. Jets that have \r]\ < 2.5 and pj > 20 GeV 
are used for the event selection, although all jets up to \r]\ < 4.5 are retained for the computation of the 
missing transverse momentum. 

Candidate stop production events, preselected by the same trigger and data quality requirements used in 
Ref. [18], are further required to contain one electron and one muon of opposite charge, with an invariant 
mass ma > 20 GeV. The leading (in pj) and next-to-leading leptons are required to have pj > 25 GeV 
and Pi > 20 GeV, respectively. 

At this stage of the selection, the background is dominated by production of top-quark pairs and Z ^ tt, 
followed by WW and Wt production. 

A requirement of mj 2 > 20 GeV, where mj 2 is the stransverse mass of the two leptons, strongly reduces 
the Z ^ TT background, which is expected to have a kinematical end-point at mT 2 = The ratio Ri 
of the and the effective mass, defined as the scalar pj sum of the the leptons and the jets, 

is useful in suppressing the tt background, which is typically characterised by a larger hadronic activity 
than in signal events. The selection chosen is > 0.3 -i- nieff (with nieff in TeV). 

After the above selections, the SM background is dominated by WW production. Two differences between 
this process and the stop pair production signal are further exploited: firstly, the WW production is dom¬ 
inated by quark-antiquark scattering, while stop pair production is mostly initiated by gluon-gluon pro¬ 
cesses, and secondly the stop pair production signal has four invisible (two neutralinos and two neutrinos) 
and two undetected (the two f^-jets) objects, while the WW process has only two. The first difference im¬ 
plies a higher longitudinal boost of the system emerging from the hard scattering in signal events than in 
background events. The variable 


AX = 2 


{Pzi^i) + Pzim 


(4) 


was defined in Ref. [119], and it is an estimator of the boost. The second difference implies a higher F"™®® 
for signal events. This is exploited by making use of 


^miss 

Ri = —^- - -. (5) 

+ /’t(^i) + Ft(^2) 

Finally, the variable cos 0b > the cosine of the angle between the direction of motion of one of the two 
leptons and the beam axis in the centre-of-mass frame of the two visible leptons [119], is sensitive to the 
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spin of the produced particles, hence it provides additional rejection power against the WW production 
process. 

A set of seven signal regions were optimised for the discovery of stop pair production, with the stop 
decaying either as tj ^ a branching ratio of 100% (assuming - m^± < 10 GeV), or as 

fj —> bW^*^x^. The definitions of the signal regions are shown in Table 4. 

Table 4; Summary of signal regions used in the analysis. The upper part of the table shows the preselection 
requirements. 


SR 

WW-SRl 

WW-SR2 

WW-SR3 

WW-SR4 

WW-SR5 

WW-SR6 

WW-SR7 

Pt(G) 

> 25 GeV 

Pt(^2) 

> 20 GeV 

Ri 

> 0.3 -t meff (TeV) 

mj2 

> 20 GeV 

AX 

<0.02 

Ri 

>0.5 

1 cos Pbl 

<0.8 

<0.8 

<0.8 

- 

- 

<0.8 

- 

nrx2 

< 45 GeV 

> 25, < 55 GeV 

- 

> 70 GeV 

> 90 GeV 

> 25, < 70 GeV 

> 80 GeV 


The background from non-prompt leptons originating from heavy-quark decays or from photon con¬ 
versions in the signal regions, or from hadrons misidentified as leptons (collectively referred to as fake 
leptons in the following), is estimated as in Ref. [18]. 

Specific confrol regions, whose evenf yield is expecfed fo be dominafed by each of fhese production 
processes, are defined and included in fhe til fo consfrain fhe normalisation paramefers. The confrol region 
CRT for tt producfion is defined by changing fhe following selections wifh respecf fo fhe signal regions: 
mj 2 > 35 GeV, Ri < 0.3. Ifs purify is 92%. The CR for WW producfion (CRW) is defined by mT 2 > 35 
GeV, AX > 0.04, and has a purify of 72%. Finally, fhe CR for Z -i- jefs (CRZ) is defined by mT 2 < 20 
GeV, 30 GeV < < 80 GeV, wifh a purify of 86%. The normalisafion facfors of fhe WW, tt, Z + jefs 

producfion processes and respectively) are defermined by a combined profile likelihood fif. 

When fesfing fhe signal-plus-background hypofhesis for rejecfion, fhe fif fakes aufomafically info accounf 
fhe signal confaminafion in fhe confrol regions. For signal scenarios considering lighf (m^- < 150 GeV) 
stops decaying fhrough t^ —> bW^*^x^^, fhe signal contamination becomes so large that jUww becomes 
unrealistically low. For such cases the fit is performed excluding CRW and taking the normalisation of 
the WW background from MC simulation. 

Systematic uncertainties, affecting both the modelling of the detector response (detector-related system¬ 
atic uncertainties) and the theoretical prediction of the cross sections and acceptances of the background 
processes (theory-related systematic uncertainties) affect the predicted rates in the signal regions. Their 
classification and estimation follows closely those defined in Ref. [18]. A few differences, discussed in 
fhe following, exisf on fhe esfimafion of fhe fheory-relafed uncerfainfies. The fofal uncerfainfy on fhe yield 
of fhe WW producfion process is composed of fhree terms: the uncertainty on the NLO hard-scattering 
calculation is taken to be the difference befween fhe predicfion of POWHEG and aMC @NLO bofh using 
PYTHIA for fhe parfon shower; fhe uncerfainfy addressing fhe choice of fhe parfon-shower model is es- 
fimafed as fhe difference of fhe aMC@NLO predicfions showered eifher wifh HERWIG or PYTHIA; fhe 
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Figure 17: Distribution of the stransverse mass mj 2 in the (a) WW-VR2 and (b) WW-VRT regions defined in the 
text. The contributions from all SM processes are shown as a histogram stack. The component labelled as “Fake 
leptons” includes the estimate of the background from non-prompt leptons. The expected signal for a model of 
stop pair production with the stop decaying as fj ^ with m,- = 160 GeV, - 150 GeV and 

m^o = 100 GeV is also shown. The lower panels show the ratio between the data and the SM prediction; the yellow 
band includes statistical and systematic uncertainties on the SM prediction. 


uncertainty due to the choice of the renormalisation and factorisation scale is evaluated by changing the 
scales independently by a factor of two or one-half and taking the maximum difference. The estimated 
relative uncertainties on the signal region yields are about 6% in SR1-SR4 and SR6; 11% in SR7 and 
29% in SR5. Similar comparisons performed on the ITZ and ZZ process yield uncertainties ranging from 
30% to 45% depending on the signal region considered. Additional systematic uncertainties are assigned 
to the small expected yields from Z -i- jets production (80%), Wt (50% to 100% depending on the SR 
considered), and non-prompt lepton background. 

The values of the normalisation factors obtained when performing the fit to the control regions only are 
shown in Table 5. 

Table 5: Normalisation factors for the ti, WW and Z -H jets background processes obtained by the combined fit to 
the control region yields. The uncertainties include systematic and statistical uncertainties. 


Normalisafion faclor 

Value 


0.94 ± 0.05 

Hww 

1.01 ±0.11 


0.95 ± 0.62 


The overall predictions of the fit are compared to the data in dedicated validation region that are kinemat¬ 
ically close to the signal region. They are defined by applying fhe preselecfion requiremenfs of Table 4 
wifh fhe addifional selecfions shown in Table 6. The mxi disfribufion in WW-VR2 and WW-VRT is 
shown in Figure 17. 

For all signal regions, fhe expecfed background yield is dominafed by producfion of WW (35% in SRI fo 
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Table 6: Summary of the validation regions used in the WW analysis. The preselection requirements of Table 4 are 
also applied in all three validation regions. 


WW-VRl 

WW-VR2 

WW-VRT 

- 

- 

0.3 < /?i < 0.3 -1- njeff ( TeV) 

0.02 < |AZ| < 0.04 

AX < 0.02 

AX < 0.02 

R2 > 0.5 

R2 < 0.5 

R2 > 0.5 

IcosPbl < 0.8 

IcosPfol < 0.8 

|cos0fo| < 0.8 



Monte Carlo (stat+syst) 
r~~l Z+jets 
□ tf 

I I ww 

Single lop (Wt) 

Fake leptons 
C3 Others 
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(a) (b) 

Figure 18: Distribution of the (a) magnitude of missing transverse momentum is™*®® and (b) stransverse mass mxa 
in WW-SR3. The contributions from all SM processes are shown as a histogram stack. The component labelled as 
“Fake leptons” includes the estimate of the background from non-prompt leptons. The expected signal for a model 
of stop parr production with the stop decaying into f, ^ bx^ —> with m,- = 160 GeV, — 150 GeV and 

m^B = 100 GeV is also shown. The lower panels show the ratio between the data and the SM prediction; the yellow 
band includes statistical and systematic uncertainties on the SM prediction. 


66% in SR4). Other important background processes are Z -i- jets in SRI (20%), non-prompt leptons in 
SR2 (12%), tt in all other SR, with contributions of about 10%. The distributions of and /nxi in the 
signal region WW-SR3 are shown in Figure 18. 

Table 7 compares the predicted and observed numbers of events in each of the signal regions. No excess 
above the SM prediction is observed, hence the results are first used to derive model-independent 95% 
CL exclusion limits on the minimum number of events beyond the Standard Model in the signal region 
assuming no signal contamination in the control regions, and then to extract limits on Cvis - cr x ex J{, 
where cr is the cross section for non-SM processes, e is the selection efficiency and ^ is fhe selection 
accepfance. These limifs are also reporfed in Table 7. Finally, 95% CL exclusion limifs are derived in 
specific supersymmefric models of direcf pair production of slops. The firsl exclusion limif (Figure 19a) 
is derived in a model where fhe slop is assumed fo decay as fj —> bx'i wilh a branching ralio of 100%, 
followed by fhe decay of fhe chargino info fhe neufralino, assumed lo be fhe sfable LSP, fhrough x\ 

chargino mass is assumed lo satisfy fhe relalion m^± = - 10 GeV, and fhe limif is 

derived in fhe nit plane. Slop masses up lo aboul 250 GeV are excluded, almosl independenfly of 

1 X 1 
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Table 7: Observed (Obs) and predicted (Exp) numbers of events in the signal regions of the WW analysis, together 
with the 95% CL upper limits on the observed and expected number of signal events and respectively), 
and on the visible cross section ((ecr)®^^). 


Signal channel 

Obs 

Exp 

c95 

‘^obs 

c9S 
^ exp 


SRI 

40 

47 ± 14 

22.6 

25.2+^;^ 

1.12 

SR2 

71 

80± 13 

25.3 

27.8+“j5 

1.24 

SR3 

215 

203 ± 27 

48.4 

46.6:« 

2.38 

SR4 

88 

81 ± 11 

35.1 

28.8+“4° 

1.73 

SR5 

4 

3.4 ± 0.9 

6.2 


0.30 

SR6 

160 

154 ± 19 

45.6 


2.25 

SR7 

21 

23 + 4 

12.4 

-3.4 

0.61 




Figure 19; Exclusion limits at 95% CL in the scenario where both pair-produced stop decay exclusively via (a) 
?! ^ followed hy x\ with lS.m{ty,x\) =10 GeV, and (b) three-body or four-body decay (depending 

on the neutralino and stop mass). The black dashed line indicates the expected limit, and the yellow band indicates 
the +lcr uncertainties, which include all uncertainties except the theoretical uncertainties in the signal. The red 
solid line indicates the observed limit, and the red dotted lines indicate the sensitivity to +lcr variations of the 
signal theoretical uncertainties. For (b), the observed limits achieved by the tlL and t2L analyses are also shown, 
and the straight dashed lines correspond to Am = m^/ -H nii, and Am = m,. 


the neutralino mass. The second limit is derived in a model where the tj decays through its three-body or 
four-body decay (depending on its mass and on that of the neutralino) into —> Mvx\ with a branching 
ratio of 100%, under the assumption that the decay happens through an off-shell top quark and an on- or 
off-shell W boson. The limit is shown in Figure 19b and fills a gap between the exclusions of the t2L and 
tlL analyses. 
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B.2.2. Final states containing two top quarks and a Higgs boson (t2tlh) 

If the lightest stop has a mass such that Am ~ nit, the sensitivity of the searches for the production 

of a tj pair is greatly reduced. One of the approaches followed is to search for direct pair production 
of t 2 instead. This is the strategy used, for example, by the t2tlZ analysis, whose signal regions were 
optimised to detect the decay of a pair-produced followed by the decay . 

Inspired by the search for a SM Higgs boson produced in association with a top quark pair, a search 
was developed and optimised for the decay > ht^, where the Higgs boson is assumed to have SM 
properties, and the tj is assumed to decay as fj —> tx\ with a BR of 100%. The final state is hence 
characterised by a large jet multiplicity, by the presence of many fj-Jets from the top quark and Higgs 
boson decays and by associated with the presence of neutrinos from semileptonic decays of the top 
quark and of neutralinos. 

The selection of electrons, muons, jets and f^-jets follows the principles outlined in Appendix A. The 
specific choices made for the px and pseudorapidity thresholds and working points of the final-state 
objects, as well as the trigger selection, are the same as those in Ref. [120]. The selection requires the 
presence of exactly one electron or muon with px > 25 GeV, F™'®* >50 GeV, at least six jets with 
px > 25 GeV and \r}\ < 2.5, of which at least two are required to be f^-tagged. The working point chosen 
for the fj-tagging is such that the efficiency to tag fj-jets (evaluated on a MC sample of tt production) is 
about 70%. 

The modelling of the production of ti pairs in association with heavy flavour (tt-i-HF) is of key relevance 
in this analysis. A detailed categorisation of tt-i-HF is made for the purpose of comparisons with different 
generators and of the propagation of systematic uncertainties on the different heavy-flavour components. 
The categorisation is also used to reweight the different flavour components of the ff-i-jets background 
to obtain a better modelling. These categorisation and reweighting procedures are discussed in detail in 
Ref. [120]. In particular, the tt + bb component, which is simulated with POWHEG, is reweighted to a full 
NLO calculation [121] performed in SHERPA 1.4.1-i-OpenLoops [122,123]. The reweighting is done at 
generator level using a number of kinematic variables such as the top quark px, tt system px, AT? and px 
of the dijet system not coming from the top-quark decay. A different reweighting is applied to the tt + cc 
and tt+ light-jets components, which is based on the ratio of the differential cross sections ai ^[s -1 TeV 
obtained in data and simulation as a function of the top quark px and tt system px [124]. 

The selected events are categorised into different channels, depending on the number of fj-tagged jets (two, 
three or at least four). The channel with at least four f^-jets has the largest signal-to-background ratio. The 
channels with two and three f7-tagged jets are used to calibrate the tt-i-jets background prediction and 
constrain the associated systematic uncertainties, which, in the channel with at least four f?-tagged jets, 
are dominated by the Z^-tagging, jet energy scale, and ff-i-jets heavy-flavour content uncertainties. 

Eor a given b-tag multiplicity, events are further categorised depending on the value of the transverse 
mass mx of the lepton and the missing transverse momentum. A “low-mx” (“high-mx”) region is defined 
by the requirement mx < 120 GeV {mj >120 GeV). 

The final discriminating variable used is defined as the scalar sum of E™®® and the transverse 

momenta of all selected jets. The signal is searched for by performing a binned likelihood fit to the 
^noiep simultaneously in the six channels defined (low/high-mx for three bins in f7-tagged jet 

multiplicity). The binning used for the distributions is that used in Eigure 20, where the background 
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Figure 20: Comparison between data and prediction for the distribution of , defined as the scalar sum of the 

missing transverse momentum and the transverse momenta of all selected jets, in the high-mx channels considered: 
(top) two h-tagged jets, (middle) three h-tagged jets and (bottom) four h-tagged jets, (left) before and (right) after 
the combined fit to data under the background-only hypothesis. The expected signal contributions from fj and 
pair production, assuming - 500 GeV, m,- = 300 GeV, = 120 GeV and a branching ratio of 100% 
for t 2 —> hfj are also shown added to the stack (red histograms, in dark red the contribution from direct fj pair 
production). The bottom panel displays the ratio of the data to the total background prediction. The hashed area 
represents the statistical and systematics uncertainty on the background. 
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estimate both before and after the fit is compared to the data in the high-mT region. The dominant post-fit 
uncertainties are those on the absolute normalisation of the tt + bb and tt + cc processes. 

The full list of detector systematic uncertainties considered, discussed in detail in Ref. [120], includes, 
beside a total uncertainty of 2.8% on the integrated luminosity, systematic uncertainties on the identifica¬ 
tion efficiency and energy scale uncertainty of the leptons, reconstruction efficiency and energy scale and 
resolution uncertainties for jets, f?-tagging efficiency and mis-tag rate uncertainties. Further modelling 
uncertainties are considered, which include, beside production cross-section uncertainties for W/Z-i-jets, 
single top and tt, dedicated uncertainties on the NLO calculation of the tt + bb process and on the model¬ 
ling of the tt + cc component. 

No significant excess above the expected background is observed, hence 95% CL limits are derived in a 
model where t 2 production is assumed, followed by the decay t^^ t^h (with a branching ratio of 100%) 
and fj —> tx\ (again with a branching fraction of 100%).^^ The limit is derived as a function of the and 
masses, under the assumption that l^m{t^,xf) = 180 GeV, and it is presented in Section 4.3. 


B.2.3. Final states containing two A-jets, a charged lepton, and missing transverse momentum (tb) 

Several phenomenological models, where both^^and^j are lighter than the stop (or the sbottom), allow 
for the tj —> tx\, t^ —> (yy* and b^ bx\, b^ —> tx\ decay channels to be open with competing branching 
ratios. Naturalness arguments require the higgsino mass parameter to be smaller than a few hundred 
GeV, while they impose virtually no constraint on the bino and wino mass parameters M\ and M 2 . If 
/r <«; Ml, M 2 , then the lightest chargino and neutralino masses are both of the order of fi and hence 
Am(x'^,X^) is small. Therefore, pair production of stops can lead to t^t^ —> tx\bx\ —> tbx^^x^^ff', where 
/ and f' represents low-px fermions emitted through—> ff'Xy Assuming both / and /' are too soft to 
be detected, the final state is characterised by the presence of a top quark, a bottom quark, and neutralinos 
escaping the detector. Similarly, b^ pair production can lead to the same final state. Dedicated SRs are 
defined that target this topology, which is not well covered by the tOL and tlL signal regions aimed at 
final states containing and the bOL signal regions targeting bbE™^^ final states. 

Both the leptonic and hadronic decays of the top quark have been studied, and the leptonic channel was 
found to give a better sensitivity to the signal models of interest. The dominant SM background processes 
in the signal regions are semileptonic tt and single top production. The SM background is evaluated using 
a combination of Monte Carlo and partially data-driven techniques. 

Events are selected online by a trigger requiring the presence of one electron or muon. The online se¬ 
lection thresholds are such that the plateau efficiency is reached for lepton transverse momenta of 25 
GeV. 

The identification criteria for electrons, muons, jets and follow the principles outlined in Ap¬ 
pendix A. In particular, electrons and muons are required to be isolated: the scalar pj sum of tracks in a 
cone AR = 0.2 around the electron (muon) is required to be smaller than 10% of the electron transverse 
momentum (1.2 GeV). The electron or muon track is excluded from the sum. The f?-tagging algorithm 
is used at an operating point with 70% efficiency in simulated top-quark pair production events. Signal 
regions are defined as detailed in Table 8, requiring one and only one electron or muon, two f?-tagged jets 
and a large Three of the SRs, labelled tb-SRin have no additional jet veto applied, while one of 

Production of q pairs is also included in the simplified models. The acceptance of the selection for such events is very small. 

Nevertheless, this component is considered as signal in the statistical analysis. 
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them (tb-SRExl) has a veto requirement on the number of jets (A^xjets) with pj >50 GeV in addition to the 
two leading Zj-tagged jets. The final SR optimisation is performed by using selections on the momenta of 
the objects, the mj and the mgff variables. In addition, the following kinematic variables are used in the 
event selection: 

- The minimum azimuthal distance between the closest b-tagged jet and the E™*®. This 
variable is used to remove multijet backgrounds with a cut of > 0.4. 

- meff! The scalar sum of the pj of the two f 7 -jets (with pj > 25 GeV and \p\ < 2.8 ) plus at most one 
light jet (with pj > 25 GeV and I 77 I < 2.5) and the The number of light jets, n, included in 
this sum depends on the signal region under study, although n = I was mostly used. 

- significance: The ratio of the to the square root of Hj, which is the scalar sum of the pj 
of the two f 7 -jets plus one light jet with px > 25 GeV and \p\ < 2.8. 

- nij: The transverse mass of the lepton and the missing transverse momentum vector. 

- rtthh: The invariant mass of the two Z?-tagged jets. 

- nihi: The invariant mass of a f 7 -tagged jet and the charged lepton. This variable is bounded from 
above at m, in tf production events. Since two jets are f?-tagged the variables nihe (1) and itihe (2) are 
defined to indicate the invariant mass constructed with the leading and subleading b-]&t respectively. 
The variable m™" is also defined to indicate the minimum between mb{ (1) and nibt (2). 

- <3/72x2: The asymmetric stransverse mass [108] is a kinematic variable which can be used to separate 
processes in which two decays giving missing transverse momentum occur. It is defined as follows: 

am\^{x)= min Uav{/72^(px(i/i),^ x^^^;T).«^T(px(t' 2 ),^x^;T)}l ( 6 ) 

where px(t^<) are reconstructed transverse momentum vectors, represent the missing transverse 
momenta from the two decays, with a total missing transverse momentum, ^x> and;^ is a free para¬ 
meter representing the unknown neutralino mass, which is assumed to be zero in the calculation. 
The 22 / 22 x 2 variable is calculated with different choices for px(yi) and px(t^ 2 )^ depending on the value 
of mbt (n) {n - 1 , 2 ), the invariant mass of the 22 * (/-tagged jet (/„ and the lepton: 

- If mbi{i) < 170 GeV and mbtij) > 170 GeV, then 22 / 22 x 2 is calculated with v\ = bi + ( and 
V2 ^ b f, 

- If mbeil) < 170 GeV and m},t{2) < 170 GeV, then 22 / 22 x 2 is evaluated using the two possible 
combinations for vi and 1 / 2 , and the minimum is used; 

- It both > 170 GeV and mbt{2) > 170 GeV the event is rejected. 

The case of both mu (1) and mu (2) exceeding 170 GeV is irrelevant: only events with the min¬ 
imum value of mbt smaller than 170 GeV populate the control, validation and signal regions. 

The optimisation is carried out using both a pMSSM signal model and simplified models where Em(x^,x\) 
5 or 10 GeV. In the case of the pMSSM model, additional non-(/-tagged jets are expected in the final state 
via the production of other SUSY particles, hence the optimisation points to SRs with no requirement on 
the Vxjets variable (tb-SRIn). In the case of the simplified models, additional jets come only from initial- 
or final-state radiation, and as a consequence a strict selection on Vxjets is applied as in the selection 
tb-SRExl. 
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Table 8: Summary of signal regions used by the tb analysis. 


SR 

tb-SRIn 1 

tb-SRIn2 tb-SRIn3 

tb-SRExl 

fj-jets 


2 fj-jets; pj >25 GeV 


1 lepton 


PT >25 GeV 




I 77 I < 2.5 (2.47) for p (e) 


(GeV) 

>200 

>120 

>220 

>160 

mj (GeV) 

>140 

>140 

>180 

>120 

meff(GeV) 

>300 

>450 

>650 

>300 

amT2(GeV) 

>180 

>200 

>180 

>180 

m™" (GeV) 


<170 


A/. 


>0.4 


significance ( GeV^'^^) 

> 8 

> 12 

>5 

> 10 

xjets 

- 

- 

<2 


The main SM backgrounds are top-pair production, W production in association with heavy-flavour jets 
and single-top production. The MC cross section is used to normalise the single-top background and all 
the other minor SM backgrounds, such as Z-i-jets, diboson production, ti+W and tt +Z. The normalisation 
factors of the tt and W + jets backgrounds are determined by a combined profile-likelihood fit. Specific 
control regions, whose event yield is expected to be dominated by each of these production processes, are 
defined and included in the fit to constrain the normalisation parameters. The tt control regions (CRT) are 
defined by inverting the selection on amT 2 , requiring amj-j <160 (180) GeV for the inclusive (exclusive) 
SRs. The purity of the tt process in the CRTs is in excess of 95%. The W + jets control regions (CRW) are 
defined by requiring mj <120 GeV. For the control regions corresponding to the tb-SRIn, events with one 
b-tagged jet are included in the CRW. Top quark pair production dominates the CRWs, with aW + jets 
purity of 30% or better. The normalisation factors luyv and are presented in Table 9. The background 
model is then validated using validation regions, where little signal contamination is expected. 


Table 9; Background scale factors for the tt and W samples, as obtained by the background fit. The errors include 
both the statistical and systematics uncertainties. 


Norm. Factor 

SRinA 

SRinB 

SRinC 

SRexA 


1.06 ±0.07 

1.12 ±0.09 

0.94 ± 0.21 

1.06 ±0.07 

Ew 

0.92 ± 0.20 

0.61 ±0.23 

0.93 ± 0.27 

1.10 ±0.34 


The distributions of the variable amj 2 in the four SRs are shown in Figure 21 together with the expected 
distribution from some of the signal models used to optimise the analysis. Table 10 compares the predicted 
and observed numbers of events in each of the signal regions. No excess above the SM prediction is 
observed, hence the results are first used to derive model-independent 95% CL exclusion limits on the 
number of events beyond the Standard Model in the signal region, and then to extract limits on cTyis = 
(T X ex tTl, where cr is the cross section for non-SM processes, e is the selection efficiency and Ji is the 
selection acceptance. All these limits are also reported in Table 10. 

Since the number of events observed agrees with the SM predictions, 95% CL exclusion limits are derived 
in specific supersymmetric models of direct pair production of stops. Simplified models were simulated 
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Figure 21: Distribution of the asymmetric stransverse mass amj 2 in the (a) SRinA, (b) SRinB (top right), (c) SRinC 
and (d) SRexA defined in the text. The contributions from all SM processes are shown as a histogram stack. The 
contribution from signal points studied by this analysis are also shown. The lower panels show the ratio between 
the data and the SM prediction; the band includes statistical and systematic uncertainties on the SM prediction. 


Table 10: Observed (Obs) and predicted (Exp) numbers of events in the signal regions of the tb analysis, together 
with the 95% CL upper limits on the observed and expected number of signal events (5®^^, and respectively), 
and on the visible cross section ((ecr)®^^). 


Signal channel 

Obs 

Exp 

c95 

obs 

c95 

exp 


SRinA 

38 

27 ±7 

28.5 

19.3!^;? 

1.41 

SRinB 

20 

14.1 ±2.8 

16.3 

10 7+^-5 

0.81 

SRinC 

10 

7.1 ±2.9 

11.9 

9.8:|3 

0.58 

SRexA 

46 

31 ±7 

32.1 

20.3!3«;° 

1.58 
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with the two decays tj —> bx^ each having a 50% BR for values of =5, 20 GeV. 

Furthermore, by using a weighted combination of these simplified models with models corresponding to 
a 100% BR in either tj —> or bx"^, limits can be obtained for any value of the stop BR. Figure 22 
shows the exclusion limits for BR(tj —> t^j)=25%, 50% and 75% for the two values of Am(xf,x^) 
considered. 

Finally, 95% CL exclusion limits are also derived for a natural pMSSM model and are presented in 
Figure 23 


C. Further details of the statistical combination of the tOL and tlL signal 
regions 

This section provides additional details on the combination of the tOL and tlL signal regions targeting 
scenarios in which the stop decays into either —> tx\ or the mixed case where tx\ and bx\ 
are both allowed, as discussed in Section 4.1 

The statistical combination of the two analyses is performed by running the combined fit simultaneously 
on the control and signal regions of the two analyses. The detector systematic uncertainties are treated as 
correlated by using, for each of the uncertainties considered, a single nuisance parameter. The supersym¬ 
metric signal parameter strength used is the same for the two analyses, while the normalisation parameters 
for the background processes are kept independent in each analysis. The nuisance parameters associated 
with modelling uncertainties of the various processes are also kept independent. 

The control regions of the two analyses are not mutually exclusive: the events that belong to both a CR 
of tOL and one of tlL are, at most about 2% of the total number of events of the tOL CR. The strategy 
adopted is to remove them from the corresponding tOL CR for the combination. It has been verified fhaf 
such removal does nol affecl fhe individual resulfs of fhe fOL analysis. 

For each combinafion performed, fhe fif sefup is validafed by checking fhaf fhe background normalisation 
parameters obfained are compafible wifh fhose obfained separately by fhe fwo analyses, by verifying fhaf 
no additional consfrainf on fhe nuisance paramefers is infroduced wifh respecf fo fhe individual fils, and 
by checking fhaf no arlificial correlation is introduced between any of the fit parameters. 

The 95% CL limit derived from the combination is shown in Figure 24, where the combined limit is 
compared to the individual limits obtained by the tOL and tlL analyses independently. 


D. Signal generation details 


Several SUSY models are considered throughout this paper. This section provides the details of how 
these signal models are generated. For all SUSY models discussed below, the detector response is sim¬ 
ulated by passing the generated events through a detector simulation [125] based on GEANT4 [126] or 
through a fast simulation using a parametric response to the showers in the electromagnetic and hadronic 
calorimeters [127] and GEANT4-based simulation elsewhere. All samples are produced with a varying 

The choice is motivated by the fact that the phase-space regions in which the two analyses determine the normalisation 
parameters of the tt, Z + jets and W + jets (for tOL) and ti and W + jets (for tlL) are characterised by different kinematic 
selections and jet multiplicities. 
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Figure 22: Exclusion limits at 95% CL from the tb signal regions for simplified models with stop decays into both 
fj —> and fj ^ bx\ and for BR(fj —» tx^)-25%, 50%, 75% (in descending rows) for the grids Am(xf,Xi) - 

5, 20 (left, right columns). The black dashed line indicates the expected limit, and the yellow band indicates the 
+ lcr uncertainties, which include all uncertainties except the theoretical uncertainties in the signal. The red solid 
line indicates the observed limit, and the red dotted lines indicate the sensitivity to +lcr variations of the signal 
theoretical uncertainties. For each point the SR giving the best expected signihcance is used. 


number of simulated minimum-bias interactions overlaid on the hard-scattering event to account for mul¬ 
tiple pp interactions in the same or nearby bunch crossings (pileup). The simulation is reweighted to 
match the number of minimum bias interactions in data, which varies between approximately 10 and 30 
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Figure 23: Exclusion limits at 95% CL from the tb signal regions for the natural pMSSM model. The black dashed 
line indicates the expected limit, and the yellow band indicates the +lcr uncertainties, which include all uncertainties 
except the theoretical uncertainties in the signal. The red solid line indicates the observed limit, and the red dotted 
lines indicate the sensitivity to +lcr variations of the signal theoretical uncertainties. For each point the SR giving 
the best expected significance is used. 


interactions in each bunch crossing. Corrections are applied to the simulated samples to account for dif¬ 
ferences between data and simulation for the trigger and reconstruction efficiencies, momentum scale and 
resolution of the final-state objects, including the efficiency of identifying jets originating from the frag¬ 
mentation of f7-quarks, together with the probability for mis-tagging light-flavour and charm quarks. 


Simplified models: The signal samples for the scenario where both stops decay to a top quark and a 
neutralino are generated using Herwig-i-i- 2.5.2 [128] interfaced to PYTHIA 6.426 [129]. The neutralino 
is fixed to be a pure bino, enhancing the decay of the component of to a right-handed top quark. 
Signal samples where the two stops decay as fj ^ bx\ are generated with MadGraph 5.1.4.8 [130]. For 
models where the W boson is on-shell, the fj decay is treated by MadGraph, while if the W is off-shell, 
PYTHIA is used fo decay fhe In fhese samples, fhe is assumed fo be mosfly a and fhe chargino 
is assumed fo decay fhrough^* —> W^*^b wifh a branching rafio of 100%. Several assumptions abouf fhe 
chargino masses are considered as described in fhe body of fhe paper. 

Models in which fhe stop is assumed fo decay eifher as tj —> fjpj or —> bx\ wifh differenf branching 

rafios are obfained by appropriafely weighfing fhree samples: one where bofh slops decay Ihrough ^ 
a second one where bofh slops decay Ihrough bx"^, and a Ihird one, where one of fhe Iwo pair- 
produced slops decays as fj ^ and fhe olher one decays as bx\- This Iasi sample is generated 
wifh MadGraph, fhe is assumed fo be a maximal mixing of and The mass of fhe chargino in fhis 
sample salisfies fhe gauge-universalify relalion m^± = 2m 

The Ihree-body slop decay samples are generaled wifh Herwig-i-i-, which performs fhe malrix elemenl 
calculation of fhe Ihree-body decay. The four-body stop decay samples are generaled wifh MadGraph. 

For all samples considered, fhe mass of fhe bolfom quark is fixed fo 4.8 GeV and ifs widlh is assumed fo 
be zero. 

The samples where bofh stops decay as cx^ are generaled wifh MadGraph, wifh one addifional 

parlon from fhe malrix elemenl. Similarly fo fhe case of fhe limif derived as a function of fhe slop 
branching rafio into tj ^ tx^, fhe samples wifh bofh stops decaying as —> cx^ and fhose where bofh 
stops decay Ihrough fhe four-body decay are appropriafely weighfed and combined wifh a Ihird sample 
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Figure 24: Combined exclusion limits at 95% CL in the scenario where both stops decay exclusively via fj ^ j. 
The black dashed line indicates the expected limit, and the yellow band indicates the +lcr uncertainties, which 
include all uncertainties except the theoretical uncertainties in the signal. The red solid line indicates the observed 
limit. For comparison the dotted green and blue lines show the expected limits from the standalone tOL and tlL 
analyses. 


where one stop decays as —> c^j and the other decays through the four-body decay to produce a sample 
of arbitrary branching ratio into —> cx^ (assuming that tj —» cx^ and the four-body decay are the only 
possible stop decays). Such mixed samples are also generated with MadGraph. 

Sbottom pair production samples are also all produced with MadGraph interfaced to PYTHIA, and no 
more than one additional parton is added to the matrix element. The PDF set used for all signal samples 
is CTEQ6L1 [131]. 


pMSSM models: In all cases, the particle spectra are generated with SOFTSUSY 3.3.3 [132], while 
sparticles decays are calculated with SUSY-HIT [133] (SDECAY 1.3b and HDECAY 3.4). The simulated 
signal events are generated using Herwig-t-i- 2.6.3 [128] with the CTEQ6E1 PDE set. 
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